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Summary 

Experimental data were obtained on an optimally contoured 
nozzle with an area ratio of 1 025 : 1 and on a truncated version of 
this nozzle with an area ratio of 440: 1 . The nozzles were tested 
with gaseous hydrogen and liquid oxygen propellants at com- 
bustion chamber pressures of 1800 to 2400 psia and mixture 
ratios of 3.89 to 6.15. This report compares the experimental 
performance, heat transfer, and boundary layer total pressure 
measurements with theoretical predictions of the current Joint 
Army, Navy, NASA, Air Force (JANNAF)developed method- 
ology. This methodology makes use of the Two-Dimensional 
Kinetics (TDK) nozzle performance code. 

Comparisons of the TDK-predicted performance to experi- 
mentally attained thrust performance indicated that both the 
vacuum thrust coefficient and the vacuum specific impulse 
values were approximately 2.0-percent higher than the turbu- 
lent prediction for the 1025:1 configurations, and approxi- 
mately 0.25-percent higher than the turbulent prediction for 
the 440:1 configuration. 

Nozzle wall temperatures were measured on the outside of a 
thin-walled heat sink nozzle during the test firings. Nozzle heat 
fluxes were calculated from the time histories of these tempera- 
tures and compared with predictions made with the TDK code. 
The heat flux values were overpredicted for all cases. The results 
range from nearly 100 percent at an area ratio of 50 to only 
approximately 3 percent at an area ratio of 975. Values ol the 
integral of the heat flux as a function of nozzle surface area were 
also calculated. Comparisons of the experiment with analyses 
of the heat flux and the heat rate per axial length also show that 
the experimental values were lower than the predicted value. 

Three boundary layer rakes mounted on the nozzle exit 
were used for boundary layer measurements. This arrangement 
allowed total pressure measurements to be obtained at 14 dif- 
ferent distances from the nozzle wall. A comparison of bound- 
ary layer total pressure profiles and analytical predictions show 


good agreement for the first 0.5 in. from the nozzle wall; but the 
further into the core flow that measurements were taken, the 
more that TDK overpredicted the boundary layer thickness. 

Introduction 

The design and analysis of efficient, high-area-ratio rocket 
nozzles requires the knowledge of core flow effects, boundary 
layer effects, contour effects, supersonic shock effects, wall 
heat transfer effects, and the specific impulse attainable. Data 
on these effects have been difficult to obtain because there are 
few altitude test facilities available for testing nozzles with area 
ratios in the range of 700: 1 to 1000: 1 . As a result, the primary 
tools for nozzle designers are theoretical methods incorporated 
in numerical codes. Many of these codes are based on the Joint 
Army, Navy, NASA, Air Force (JANN AF) prediction method- 
ology outlined in reference 1. One of the computer programs 
most often used for nozzle analysis is the Two-Dimensional 
Kinetics (TDK) Nozzle Performance Computer Program 
(ref. 2). As stated in reference 3, when the methodology was 
developed, area ratios of 1 00: 1 were considered large-area ratio 
nozzles. In the past 20 years, the need for increased perfor- 
mance from orbital transfer vehicles has required investigation 
into area ratios of up to 1 000: 1 . Several experimental programs 
have been undertaken to validate codes developed with the 
JANNAF methodology at higher area ratios (refs. 3 to 7). As a 
result of these activities, the codes are considered validated for 
low-area-ratio nozzles (up to 300:1) and are being used to 
extrapolate results to high-area-ratio nozzles. These extrapola- 
tions lack confidence without experimental validation and raise 
questions as to the relevance of trade studies for iuture rocket 
engine designs. Hence an experimental program (ref. 8) was 
undertaken to provide data to validate the codes for high-area- 
ratio nozzles at high chamber pressures. As part of this effort, 
a series of tests were conducted in the altitude test capsule at the 
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Figure 1 . — NASA Lewis Research Center’s Rocket Engine Test Facility. 
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NASA Glenn Rocket Engine Test Facility (RETF). Previous 
tests in this program were reported to be in the laminar houndary- 
layer regime (refs. 3 to 5) at a nominal combustion-chamber 
pressure of 2.4 MPa (350 psia) and at a Reynolds number, based 
oh throat diameter, that ranged from 3.1 IxlO 6 to 4.14x10 s . 
Present tests were considered to be in the turbulent boundary- 
layer regime at combustion-chamber pressures that ranged 
from 12.4 to 16.5 MPa (1800 to 2400 psia) and Reynolds 
numbers, based on throat diameter, that ranged from 1 43x 1 0 6 
to 2.05x 1 0 6 (ref. 9). The nozzles used in these tests had nominal 
2.54-cm- ( 1 .00-in.- ) diameter throats with area ratios of 1 025: 1 
and 440:1, and were fired with gaseous hydrogen and liquid 
oxygen. This report compares the performance and heat trans- 
fer test results with the theoretical predictions of the TDK 
computer code. In addition, boundary layer rakes were used to 
measure the total pressure profile of the boundary layer for 
comparison with analytical predictions. A symbols list is pro- 
vided in appendix A. 

Apparatus 

Facility 

Testing was conducted at the NASA Glenn Rocket Engine 
Test Facility (RETF) (fig. 1) and utilized on the facility’s 


altitude test capsule, thrust stand, propellant feed system, and 
data acquisition system. The altitude test capsule (fig. 2) 
simulated the static pressure at altitude by three methods of 
vacuum pumping, all acting simultaneously. The first method, 
a second-throat diffuser, utilized the kinetic energy of the 
rocket exhaust to pump the nozzle flow into a spray cooler. The 
second method chilled the exhaust gas in the spray cooler where 
approximately half was condensed to liquid water and drained. 
The third method pumped the remaining uncondensed exhaust 
by nitrogen-driven ejectors. The facility ejector system reduced 
the capsule pressure to approximately 4. 1 kPa (0.6 psia), with 
further pumping accomplished by the engine exhaust. Addi- 
tional facility details are given in references 4 and 8. 

The thrust stand, which had a full-scale measurement range 
of 1 7.8 kN (4000 lb f ), w as designed to have a standard deviation 
(2-o) variation of less than ±0. 1 percent of full scale. With the 
test capsule at altitude pressure, the thrust stand w'as calibrated 
against a reference load cell, which had a 2-0 variation of less 
than ±0.05 percent of full scale and a calibration traceable to 
the National Institute of Standards and Technology (NIST). 

The propellant feed system consisted of a gaseous hydrogen 
fuel circuit and a liquid oxygen oxidizer circuit. High-pressure 
gaseous hydrogen bottles comprised the fuel circuit; the oxi- 
dizer circuit was a high-pressure liquid oxygen tank pressur- 
ized from high-pressure gaseous helium bottles (fig. 3). The 
flow' rates were measured with calibrated venturis. 


Exhaust 
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D1 Fuel venturi AP 
D2 Oxygen venturi AP 
D3 Fuel injection A P 
D4 Oxygen injection A P 
D5 Nozzle wall A P 
D6 Nozzle wall A P 
D7 Nozzle wall AP 
D8 Altitude AP 


D9 Altitude AP 
Ft Thrust 

PI Fuel supply pressure 
P2 Oxygen supply pressure 
P3 Fuel injection pressure 
P4 Oxygen injection pressure 
P5 Chamber pressure 
P6 Chamber pressure 


P7 Vacuum reference pressure 
T1 Fuel supply temperature 
T2 Oxygen supply temperature 
T3 Fuel injection temperature 
T4 Oxygen injection temperature 
T5 Nozzle wall temperature 
T6 Nozzle waJI temperature 
T7 Nozzle wall temperature 


Figure 3. — Schematic showing propellant circuit and instrumentation. 



Figure 4. — Rocket engine injector. 


Injector- 


Radius, 
2.54 (1) 



/ 


\ 


Combustion chamber 
pressure tap 


Figure 5. — Schematic of rocket combustion chamber. 
All dimensions in inches (centimeters). 


Test Hardware 

The test hardware consisted of an injector, chamber, nozzles, 
and boundary layer total pressure rakes. The injector (fig. 4) had 
a porous faceplate for gaseous hydrogen injection and 36 tubes 
for liquid oxygen injection. A gaseous hydrogen and gaseous 


oxygen torch igniter located in the center of the injector ignited 
the propellant mixture. As shown in figure 5, the combustion 
chamber was a water-cooled copper spool 15.24-cm (6- in.) 
long with an inside diameter of 5.22 cm (2.055 in.). 

Two low-area-ratio nozzles, e = 1 0.7: 1 and 4: 1 (fig. 6), were 
used to calibrate the effective combustion chamber pressure 
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Figure 6. — Sea-level engine installed on test stand. 



Figure 7. — High -area-ratio nozzle on test stand. 


Axial distance 
from throat 

Radius 

cm 

in. 

cm 

in. 

0 

0 

1 .2700 

0.5000 

.3929 

.1547 

1 .4371 

.5658 

.4641 

.1827 

1 .4961 

.5890 

.6068 

.2389 

1.6190 

.6374 

.7503 

.2954 

1 .7404 

.6852 

.8230 

.3240 

1.8031 

.7099 

1.3246 

.5215 

2.2426 

.8829 

1.7844 

.7025 

2.6515 

1 .0438 

2.3777 

.9361 

3.1643 

1 .2458 

3.2062 

1 .2623 

3.8572 

1.5186 

7.0256 

2.7660 

6.6703 

2.6261 

7.8931 

3.1075 

7.2426 

2.8514 

9.6269 

3.7901 

8.3320 

3.2803 

10.6505 

4.1931 

8.9433 

3.5210 

11.6738 

4.5960 

9.5341 

3.7536 

12.9022 

5.0796 

10.2189 

4.0232 

15.3429 

6.0405 

11.5108 

4.5318 

16.5392 

6.5115 

12.1150 

4.7697 

19.5651 

7.7028 

13.5702 

5.3426 

23.3688 

9.2003 

15.2710 

6.0122 

25.4869 

10.0342 

16.1651 

6.3642 

29.5410 

11.6303 

1 7.7871 

7.0028 

33.7297 

13.2794 

19.3558 

7.6204 

36.2996 

14.2912 

20.2705 

7.9805 

38.8696 

15.3030 

21 .1 524 

8.3277 

41.4193 

16.3068 

21 .9977 

8.6605 

47.2194 

18.5903 

23.8201 

9.3780 

51.1703 

20.1458 

24.9895 

9.8384 

55.1213 

21 .7013 

26.1064 

10.2781 

60.4944 

23.8167 

27.5486 

10.8459 

71.1091 

27.9957 

30.1694 

11.8777 

76.2211 

30.0083 

31 .3365 

12.3372 

90.6396 

35.6849 

34.3444 

13.5214 

105.3071 

41.3532 

36.9933 

14.5643 

113.0838 

44.5212 

38.3365 

15.0931 

128.5725 

50.6191 

40.6598 

16.0078 


Carbon 



Figure 8. — Thruster assembly showing components and 
expansion area ratios, e. 
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Axial distance from throat, in. 

Figure 9. — Nozzle contour and coordinates. 
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Location 

Nominal 
distance from 
nozzle wall, 
in. 

Actual 

distance from 
nozzle wall, 
in. 

1037:1 Area ratio nozzle 

A 

0.0655 

0.0655 

B 

.6855 

.6795 

C 

1 .3055 

1.3135 

D 

1 .9255 

1.9345 

E 

2.5455 

2.5645 

F 

3.1655 

3.1615 

G 

3.7855 

1 

3.7865 

440:1 Area ratio nozzle 

A 

0.0900 

0.0900 

B 

.7100 

.7045 

C 

1.3300 

1 .2970 

D 

1.9500 

1.9155 

E 

2.5700 

2.5375 

F 

3.1900 

3.1515 

G 

3.8100 

3.7625 


(a) 



Location 

Nominal 

Actual 


distance from 

distance from 


nozzle wall, 

nozzle wall, 


in. 

in. 

1037:1 Area ratio nozzle 

A 

0.0440 

0.0440 

B 

--- 

— - 

0 

.4040 

.3845 

D 

--- 

— 

E 

c 

.7640 

.7535 

r 

G 

1.1240 

1.1280 

440:1 Area ratio nozzle 

A 

0.0680 

0.0680 

B 

— 

— 

C 

.4280 

.4195 

D 

— 

--- 

E 

.7880 

.7455 

r 

G 

1.1480 

1,0805 


Figure 10. — Boundary layer rakes, (a) Nominally 4-in. high, (b) Nominally 1-in. high. 
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Location 

Nominal 
distance from 
nozzle wall, 
in. 

Actual 

distance from 
nozzle wall, 
in. 

1037:1 Area ratio nozzle 

A 



— 

B 

0.0895 

0.0895 

C 

— 

— 

D 

.4495 

.4380 

E 

F 

.8095 

.8240 

G 

— 

— 

440:1 Area ratio nozzle 

A 





B 

0.2845 

0.2845 

C 

--- 

— 

D 

.6445 

.6440 

E 

... 

— 

F 

1.0045 

1.0125 

G 

— 

— 


(C) 


Figure 10.— Concluded, (c) Alternate nominally 1 -in. high. 


P at the nozzle entrance as a function of the static pressure 
P at the end of the combustion chamber. Two high-area-ratio 

CM _ 

nozzle configurations, £ = 1025: 1 (fig. 7) and 440: 1, were used 
to obtain performance data. The nozzle converging-diverging 
section was a water-cooled copper throat piece that started at 
the 5.22-cm (2.055-in.) combustion chamber inside diameter, 
converged to the 2.54-cm (1.0-in.) throat, and diverged to an 
area ratio of 30: 1 . At this point, a 0.635-cm- (0.25-in.-) thick 
carbon-steel nozzle skirt was attached that continued the con- 
tour to an expansion area ratio of 440: 1 . The final piece of the 
nozzle was a 0.635-cm- (0.25-in.-) thick carbon-steel skirt 
extension that concluded the contour to an area ratio of 1025: 1 
(fig. 8). The carbon-steel sections of the nozzle skirt, which 
were not actively cooled, were designed to survive the exhaust 
temperatures by nature of their inherent heat sink. 

Contour calculations were based on the Rao nozzle optimi- 
zation process (ref. 10). which uses a Rao nozzle design code 
(ref. 1 1 ). Figure 9 shows a plot and a table of the nozzle 
coordinates. A row of static pressure taps through the wall of the 
carbon-steel nozzle skirt measured the nozzle wall static pres- 
sures, and chromel-constantan thermocouples spot-welded to 
the outside surface measured the temperature of the outside 
wall of the carbon-steel skirts. These thermocouples were pre- 
referenced to a 67 °C (150 °F) oven. Their specified absolute 
accuracy was ±1.1 K (±2 °F). Temperatures were measured 
at nine axial locations in a row, circumferentially displaced 45° 
from the static pressure tap locations. 



Figure 1 1 . — Nominal 1 -in. rake mounted on a nozzle. 


Boundary-layer total pressure measurements were made 
with a series of three total pressure rakes (figs. 1 0(a) to (c)) that 
were constructed with massive copper bodies to provide con- 
duction and a heat-sink for the main probe support structure. 
The individual probe tubes were made of 0.203-cm- (0.08-in.-) 
diameter tubing to provide adequate spatial resolution to the 
pressure profiles, yet they were not so fine as to have no thermal 
survivability. These tubes were made of a moly-rhenium alloy 
to provide some additional thermal survival capability. Figure 1 1 
shows a nominal 2.54-cm (1-in.) rake mounted on a nozzle. 
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Procedure 

Experimental Procedure 

Atmospheric Testing . — Atmospheric pressure tests were 
first performed with the two low-area- ratio nozzles (e = 10.7: 1 
and 4:1) to determine P . The firings were approximately 
3 sec in duration. A steady-state condition was reached at or 
before 2.5 sec; this provided about 0.5 sec of steady-state 
operation before shutdown. 

Altitude Testing . — The high-area-ratio nozzles (£ = 1025: 1 
and 440:1) were tested at altitude. A typical attitude firing 
started with the gaseous nitrogen ejectors evacuating the test 
capsule and spray cooler to a pressure of approximately 4. 1 kPa 
(0.6 psia). At this pressure, the thruster was fired for about 3 sec. 
The pumping action during firing further reduced the pressure 
in the test capsule from 4.1 to approximately 1.4 kPa (0.6 to 
~0.2 psia). A steady-state pressure condition was reached at, 
or before, 2.5 sec, again providing about 0.5 sec of steady-state 
operation before shutdown. 

At thruster shutdown, the exhaust How through the diffuser 
stopped, and a pressure pulse propagated from the spray cooler 
to the test capsule, raising its pressure to the original 4.1 kPa 
(0.6 psia). Simultaneously, the two isolation valves between 
the ejectors and the spray cooler were closed and the ejectors 
were turned off. 

The high-area-ratio nozzles (1037:1 and 440:1) with the 
boundary layer rakes were tested at altitude exactly as the high- 
area-ratio ( 1 025: 1 and 440: 1 ) nozzles without rakes, except that 
some operational techniques were added to improve the surviv- 
ability of the rakes. The addition of a new throat section, which 
was required for testing at the higher area ratio with boundary 
layer rakes, resulted in the 1037:1 configuration. Conventional 
transducer installation at the end of some tubing length from the 
rake would require the tubing and transducer volume to till with 
combustion gas until the pressure reached equilibrium. This 
How into the rake would have transferred significant heat to the 
thin wall sections of the probe and tubing, resulting in melting 
or burning. Such inflow was avoided by filling the transducer, 
the connecting tubing, the rake, and the probe with room 
temperature gaseous nitrogen at a pressure of very nearly full 
scale on the transducer. This pressure produced a continuous 
outflow through the rake and probe tube that achieved three 
beneficial effects: (1) hot gas would not flow into the rake, 
(2) all the attendant hardware was cooled convectively by the 
outflow, and (3) the outside of the rake was shielded and/or 
film-cooled by spillage of the out-flow over the outside of the 
rake. No boundary-layer measurements could be made during 
the gaseous nitrogen purge. However, the rake was well pro- 
tected during thruster startup and until the nozzle flow achieved 
steady-state conditions. Once at steady-state, the purge flow was 
stopped by a high-speed solenoid valve. Then, the gaseous 
nitrogen bled down until it w as at the same pressure as the nozzle 
total pressure at the tip of the rake probe. This was then 


recorded, and the thruster was shut down. The entire duration of 
the shutoff gaseous nitrogen rake purge was 0.5 sec. This was 
sufficient to allow the transducer to get well into steady-stale 
pressure. 

Analytical Procedure 

Experimental results for all the tests were compared with 
analytical predictions from the Liquid Propulsion Program 
(LPP) (June 1994) version of the TDK code. This program 
performs two-dimensional equilibrium, frozen, or kinetic nozzle 
performance calculations with boundary layer effects (ref. 2). 
The computational portion of TDK consists of six modules: 
one-dimensional equilibrium (ODE), one-dimensional kinetic 
(ODK), transonic flow (TRANS), method of characteristics 
(MOC), and two boundary layer modules (BLM and MABL). 
Figures 1 2(a) and (b) show the distribution of the modules in the 
nozzle along with a master flowchart of the program (ref. 12). 
A brief description of the modules follows. Additional infor- 
mation can be found in references 2, 3, and 12. 

The ODE module calculates one-dimensional ideal rocket 
engine performance using either chemical frozen or chemical 
equilibrium conditions. The ODK module calculates in viscid 
one-dimensional equilibrium, frozen, and nonequilibrium nozzle 
expansion of gaseous propellant exhaust flows. The TRANS 
flow module calculates two-dimensional flow conditions in the 
transonic region of the nozzle throat. TDK uses this informa- 
tion toobtainan initial data line for the MOC module. The MOC 
module calculates the loss in nozzle performance caused by 
flow divergence. A finite difference mesh was constructed by 
tracing gas streamlines and characteristic surfaces. A separate 
boundary layer analysis was performed by using both the BLM 
and the MABL modules. As reported previously (ref. 5), the 
BLM module calculates compressible laminar and turbulent 
wall boundary layers in axisymmetric nozzles. BLM uses the 
Keller and Cebeci (ref. 13) two-point finite difference method 
to calculate the boundary layer properties and uses the Cebeci- 
Smith (ref. 1 4) eddy-viscosity formulation to model the turbu- 
lent boundary layer. 

The MABL module found in TDK is a modified version of 
the original MABL module, which was developed in 1971 by 
Levine (ref. 15). Unlike BLM, MABL allows users to chose 
either shifting equilibrium, frozen chemistry, or finite rate 
kinetics to govern the boundary layer flow chemistry. In the 
current analysis, the code was run with finite rate kinetics for 
the MABL module. As with BLM, the Cebeci-Smith eddy- 
viscosity model is used to model the turbulent boundary layer. 

Both hardware dimensions and experimental test conditions 
were input to the TDK code to model nozzle performance. 
Table I gives the geometry input for the combustion chamber 
section through the tangent point of the throat exit radius, and 
figure 9 gives the coordinates for the nozzle contour. Each 
point was normalized by the throat radius before being input to 
the program. Table II shows the experimental values input to 
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Throat 


MABL 



kinetic 

( a ) analysis 


Main program 



i 
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Figure 12. — TDK analysis (ref. 2). (a) Schematic, (b) Master flow chart. 


TAHl.K I. — UHOMKTRY INPUT TO TDK H)R C'( JMIHJSTK )N 
CHAMBER SECTION 


Parameter 

TDK variable 

fix pansit 

>n area ratit>. t' 

1025:01 

440: 1 

Throat radius, cm (in.) 

RSI 

1.27 (0.5) 

1.262 (0.497) 

Inlet contraction ratio 

HC RAT 

4.22} 

4.274 

Inlet wall radius’' 

RI 

2 

2 

Inlet angle. deg 

THHTAI 

25 

25 

Upstream wall radius oi curvature 1 

KWTlf 

2 

2 

Downstream wall radius of curvature' 1 

RWTI ) 

0.4 

0.4 

No/./lo attachment angle, deg 

THI-TA 

.59.4 1 

59.41 

No/.zlo exit angle, deg 

THh 

7.94 

15.5 


‘Normalized by throat radius. 
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TABLH II — RESULTS OF ALTITUDE PRESSURE THSTS 


Reading 

Nozzle throat 
area. A t 

Nozzle exit 
expansion 
area ratio, 

k' 

Measured chamber pressure 

Effective chamber 
pressure, a 
p c.c 

Propellant 

mixture 

ratio, 

O/F 


ei»“ 

in." 



At injector end. 

Corrected for 











P.V., 

momentum 













pressure loss. 













P 

c,T 











MPa 

psia 

MPa 

psia 

MPa 

psia 


569 

5.067 

0.7854 

1025 

12.485 

1810.8 

12.448 

1805.3 

12.326 

1787.7 

3.89 

570 







12.867 

1866.1 

12.797 

1856.0 

12.645 

1834.0 

5.97 

571 







12.675 

1838.3 

12.621 

1830.4 

12.488 

1811. 1 

4.70 

575 







14.562 

2111.9 

14.502 

2103.3 

14.350 

2081.2 

4.65 

576 







14.850 

2153.8 

14.775 

2142.9 

14.605 

2118.2 

5.68 

577 

1 


1 


1 


14.429 

2092.7 

14.373 

2084.6 

14.225 

2063.1 

4.47 

580 







16.586 

2405.5 

16.531 

2397.5 

16.364 

2373.3 

4.27 

601 

5.007 

.7760 

440 

12.993 

1884.4 

12.923 

1874.3 

12.768 

1851.8 

6.15 

602 

5.(X)7 

.7760 

440 

12.740 

1 847.7 

12.681 

1839.2 

12.542 

1819.0 

5.11 

603 

5.1X17 

.7760 

440 

12.621 

1830.4 

12.581 

1824.7 

12.457 

1806.7 

4.01 


Calculated with low nozzle exit expansion area ratio e correlation. 


Reading 

Vacuum thrust, 

p'v 

Ambient 
pressure around 
nozzle, 

P a 

Characteristic 

exhaust 

velocity, 

C* - 

Characteristic 
exhaust velocity 
efficiency, 

1c » 

percent 

N 

tor 

kPa 

psia 

m/s 

ft/s 

569 

1 1 863 

2667.1 

1.491 

0.2162 

2476 

8124 

98.9 

570 

12 957 

2913.0 

1.342 

.1947 

2330 

7643 

98.6 

571 

12 392 

2785.9 

1.313 

.1905 

2448 

8033 

99.7 

575 

14 179 

3187.7 

1.470 

.2132 

2448 

8033 

99.5 

576 

14 904 

3350.8 

1.510 

.2190 

2372 

7782 

99.4 

577 

14010 

3149.8 

1.446 

.2097 

2467 

8094 

99.8 

580 

16 109 

3621.7 

1.582 

.2295 

2490 

8170 

1(X).2 

601 

12 498 

2803.7 

.9143 

.1326 

2328 

7637 

99.2 

602 

1 1 923 

2680.5 

.7812 

.11 33 

2416 

7925 

99.5 

603 

1 1 450 

2574.1 

.6943 

.1007 

2497 

8192 

100.0 


Reading 

Fuel injection 

Oxidizer injection 

Propellant 
flow' rate. 

rii 

Pressure, 

Temperature, 

T fi 

Pressure. 

P oi 

Temperature. 

T 

01 


M Pa 

psia 

K 

°R 

MPa 

psia 

K 

°R 

kg/s 

Ibm/s 

569 

16.563 

2402.2 

297 1 

534.8 

13.509 

1959.3 

112.6 

202.6 

2S22 

5.561 

570 

15.316 

2221.3 

297.1 

534.8 

14.393 

2087.4 

117.8 

212.1 

2.751 

6.064 

571 

15.863 

2300.7 

297.3 

535. 1 

13.967 

2025.6 

121.6 

218.8 

2.584 

5.697 

575 

18.317 

2656.6 

296.3 

533.3 

16.138 

2340.6 

108.6 

1954 

2.970 

6.547 

576 

17.837 

2586.9 

296.8 

534.2 

16.778 

2433.3 

111.6 

200.9 

3.120 

6.878 

577 

18.353 

2661.8 

296.8 

534.3 

15.998 

2320.3 

115.0 

207.0 

2.922 

6.441 

580 

21.422 

3106.9 

298.9 

538. 1 

18.521 

2686.1 

106.6 

191.8 

3.329 

7.340 

601 

15.31 1 

2220.6 

3(X).7 

541.3 

14.480 

2100. 1 

109.1 

196.3 

2.746 

6.054 

602 

15.570 

2258.2 

299.5 

539. 1 

14.011 

2032.1 

109.6 

197.2 

2.600 

5.731 

603 

16.431 

2383.1 

299.3 

538.8 

13.707 

1987.9 

113.0 

203.4 

2.498 

5.506 


Reading 

Measured 
vacuum thrust 
coefficient, 

C KV 

Thrust 

coefficient 

efficiency, 

"cp.v 

percent 

Vacuum 

specific 

impulse. 

f sp,V- 

s 

Vacuum specific 
impulse efficiency. 

^sp.V ’ 
percent 

569 

1.900 

96.8 

479.6 

95.8 

570 

2.022 

96.3 

480.4 

95.0 

571 

1.958 

97.3 

489.0 

96.9 

575 

1.950 

97.1 

486.9 

96.9 

576 

2.014 

97.0 

487.2 

96.4 

577 

1.944 

97.3 

489.0 

97.1 

580 

1.943 

97.9 

493.4 

98.2 

601 

1.955 

94,0 

464.1 

93.2 

602 

1.899 

94.2 

467.7 

93.7 

603 

1.836 

94.2 

467.5 

94.1 


)() 
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TABLE 111 — NOZZLE INSIDE WALL TEMPERATURES 


Reading 

Effective combustion 
chamber total pressure 
at nozzle entrance. 

fVe 

Propellant 

mixture 

ratio. 

O/F 

Expansion area ratio, e 

so - 

50.6 

100 1 

Nozzle wall temperature | 

MPa 

psia 

K 

°R 

K 

°R 

K 

°R 

569 

12.326 

1787.7 

3.89 

431.55 

776.79 



363.59 

654.47 

570 

12.645 

1834.0 

5.97 

518.48 

933.26 



419.70 

755.46 

571 

12.488 

181 1.1 

4.7 

491.77 

885.18 



407.41 

733.33 

575 

14.350 

2081.2 

4.65 

503.29 

905.92 



405.86 

730.55 

576 

14.605 

2118.2 

5.68 

536.22 

965.19 


— 

434.23 

781.62 

577 

14.225 

2063.1 

4.47 

495.39 

891.71 



421.26 

758.27 

580 

16.364 

2373.3 

4.27 

520.37 

936.67 



418.81 

753.85 

601 

12.768 

1851.8 

6.15 



505.34 

909.62 


— 

602 

12.542 

1819.0 

5.1 1 



496.83 

894.30 



603 

12.457 

1806.7 

4.01 



470,64 

847.16 




Reading 


101.2 


K 


569 


570 


571 


575 


576 


577 


580 


601 

■ 

602 


603 



Expansion area ratio, e 
2(X) | 202.4 | 300 

Nozzle wall temperature 


°R 


K 

°R 

326.67 

588.01 

365.39 

657.70 

360.8 1 

649.45 

350.89 

631.61 

378.03 

680.46 

378.01 

680.41 

359.93 

647.88 



R 


K 

314.03 

345.37 

344.37 
330.71 

356.04 
362.50 
337.52 


°R 

565.25 

621.67 

619.87 
595.28 

640.88 
652.50 
607.54 


K 


303.6 


°R 


744.41 

782.41 
755.79 


354.56 

387.57 

386.30 


638.21 

697.63 

695.34 


335.68 

367.03 

373.26 


604.23 

660.66 

671.87 


Reading 

Expansion area ratio, e. 


388.0 1 

392.7 

500 

635 

800 

975 1 


Nozzle wall temperature | 


K 

°R 

K 

°R 

K 

°R 

K 

°R 

K 

°R 

K 

°R 

569 

312.41 

562.34 



306.48 

551.66 

309.00 

556.20 

310.91 

559.63 

308.89 

556.01 

570 

336.46 

605.63 



326.89 

588.40 

322.96 

581.32 

320.07 

576.13 

314.78 

566.60 

571 

337.66 

607.78 



329.03 

592.25 

327.1 1 

588.79 

326.74 

588.14 

323.06 

581.51 

575 

324.82 

584.67 



314.52 

566.14 

313.56 

564.40 

311.21 

560.17 

310.17 

558.31 

576 

347.07 

624.73 



335.72 

604.30 

329.56 

593.20 

325.63 

586.14 

320.85 

577.53 

577 

356.27 

641.28 



344.16 

619.48 

340.23 

612.42 

337.56 

607.61 

334.11 

601.39 

580 

329.67 

593.41 



319.46 

575.02 

317.65 

571.77 

314.06 

565.31 

313.92 

565.06 

601 



331.71 

597.08 









602 



360.11 

648.20 









603 



370.24 

666.44 










TDK: effective chamber pressure, mixture ratio, fuel injection 
temperature, and oxidizer injection temperature. The propel- 
lant injection temperatures were used to determine the propel- 
lant enthalpies. Since the fuel was gaseous hydrogen, the TDK 
code used the enthalpy that was based on the input tempera- 
ture. The GASPLUS code (ref. 16) was used to determine the 
enthalpy for liquid oxygen on the basis of the experimentally 
determined inlet pressure and temperature to the injector. 
Because GASPLUS has a different reference state than TDK. 
the enthalpy values were corrected for this. 

The experimentally determined outside wall temperatures 
and their time rate of change were used to calculate the inside 


nozzle wall temperatures according to the method described in 
reference 5. These calculated temperatures were then used as 
nozzle input (table III). Conditions on the water-cooled com- 
bustion chamber wall from the injector face through the throat 
plane were not available and had to be estimated. For the purpose 
of this analysis, wall temperatures were distributed between 
700 and 844.4 K (1260 and 1520 °R) in this region. This 
temperature range was selected on the basis of previous com- 
bustion chamber testing data (refs. 17 and 18). Sensitivity ot the 
TDK code results to variations in combustor wall temperatures 
was negligible, especially in comparison to the effects of 
chamber pressure and mixture ratio on the final results. 
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The TDK code requires that the boundary layer be set to 
either laminar, turbulent, or transitional How at either a specific 
nozzle location or at a specific Reynolds number based on the 
momentum thickness Re$. When Re§ = 400, the boundary layer 
transitions to turbulent How (ref. 5). Although it was assumed 
that the boundary layer was turbulent, both TDK/BLM and 
TDK/MABL were run in transition with Re$ = 400 to allow the 
code to estimate the exact transition point. In all cases, the code 
indicated that transition occurred near the injector face in the 
combustion chamber; therefore, the program results were based 
on a fully turbulent boundary layer. Asapointofcomparison, but 
of no physical significance, the code was also run with a lami- 
nar boundary layer assumption. These results are reported in 
appendix B. 

The analytically predicted total pressures in the rakes were 
obtained from the LPF version of the TDK code using a pitot 
probe subroutine. This subroutine provides total flow condi- 
tions to simulate the placement of a pitot probe into the flow 
field at specified radial locations. 


Experimental Data Analysis 

Performance 


where F is the aneroid-corrected thrust and A is the nozzle 
exit area. 

Effective Chamber Pressure . — For the effective combus- 
tion chamber total pressure at the nozzle entrance P to be truly 
representative, a thorough survey of the distribution of pres- 
sures in the combustion chamber would have had to have been 
made by taking readings from several static pressure laps in the 
combustion chamber. Then, these measurements would have had 
to have been integrated and averaged to obtain an integrated 
mean pressure that could be corrected for momentum pressure 
loss and used as P . In an alternative method that was used for 
the present study, P was determined by the following equation: 


1 

P ( T ' 

' P i.e ' 

II 

or 

k.J 

UrJ 


(3) 


where P is the chamber pressure measured at a single in jector 
faceplate position, P. y7P. f/ is the conversion of the chamber 
static pressure before combustion to total pressure after com- 
bustion (momentum pressure loss), and P IP ^is the correc- 
tion that accounts for any variations in pressure distribution 
across the injector face. The momentum pressure loss was 
calculated by the following equation from reference 20: 


Propellant Mass Flow . — Propellant mass flows were meas- 
ured with calibrated venturis. Each mass flow was calculated 
from conditions at the venturi throat by 



' P* | ftr - Kn- " 

y P T c e < , 


(4) 


m-C tt pA v V (I) 

where C {! is the venturi discharge coefficient, p is the throat 
density, A is the venturi throat area, and V is the velocity; p 
and V were calculated from one-dimensional mass and energy 
equations; and real fluid properties were obtained from the 
fluid properties program GASP (ref. 19). Venturi calibrations 
of Cy were performed by the Colorado Engineering Experi- 
ment Station. Values of the discharge coefficient were trace- 
able NIST, and the uncertainty values were ±0.5 percent 
of full scale. 

Vacuum Thrust . — The vacuum thrust was determined by 
measuring the thrust produced at the test capsule ambient pres- 
sure P and by applying two corrections. The first correction 
compensated for the thrust-stand zero shift that occurred from 
the change in capsule pressure during thruster startup. This 
correction, referred to as an aneroid correction, is explained in 
reference 4. The second correction adjusted the thrust meas- 
ured at a P a of approximately 1 .4 kPa (0.2 psia) to the thrust that 
would have been measured if P a had been an absolute vacuum. 
This thrust was calculated by adding the force induced by the 
capsule pressure on the nozzle exit area to the measured thrust: 

F v = F + {P a xA ex ) (2) 


where P/P r is the static-to-total pressure ratio in the combus- 
tion chamber; / is the theoretical subsonic specific impulse 
inside the combustion chamber; g is the proportionality con- 
stant; V uv is the propellant mass-averaged injection velocity; 
C*Tii(()l)l:) * s theoretical characteristic exhaust velocity, and 
£ r is the thruster contraction area ratio. The ratio P IP T was 
derived semi-empirically by the following procedure. A series 
of low-area-ratio nozzle tests (£= 10:7 and 4: 1 ) were performed 
to develop a correlation between single-point chamber pressure 
measurements corrected for momentum pressure loss and the 
effective chamber pressure. These two pressures are defined at 
the same axial location in the chamber and vary only in that P T 
defines a single point and P.. e defines an average pressure at 
that axial location. This procedure is a calibration of the injector 
and chamber pressure tap. In these tests, the contour of the 
combustion chamber up to the throat was identical to that used 
in the test of the high-area-ratio nozzles. 

The contour downstream of the throat was identical to that of 
a low-area-ratio divergent nozzle with a thrust coefficient 
calculated by an iterative procedure using the TDK program. 
The calculated thrust coefficient obtained from TDK was used 
with the experimental measurements of thrust from the low- 
area-ratio tests and with the P.. e calculated by the following 
equation: 


12 
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^F.V,ThiTDK) A t 

where C FVTh[TDK) is the theoretical, two-dimensional- 
kinetics, vacuum thrust coefficient and A f is the nozzle throat 
area. Next, the values of P. e were related to the calculated total 
pressure after combustion P c f , and a correlation was devel- 
oped. This correlation, P ( e IP. T , which was plotted versus the 
propellant mixture ratio OIF , represents the correction lor 
nonuniform pressure distributions (fig. 1 3). A straight line was 
fit to the data with a least-squares best fit, and the equation of 
this line was used as the correlation. 

Equation (3) is valid because the same injector and chamber 
contour were used in both the low-area-ratio and high-area- 
ratio tests. The chamber static pressure was measured at the 
injector face static tap to obtain P C(f , and the momentum 
pressure loss conversion (eq. (4)) provided P c j!P C (jr The semi- 
empirical correlation P c JP c j versus OIF Irom the low-area- 
ratio nozzle tests provided the P { e /P c ^correlation. 

Performance Calculations . — By definition. 


The values of P , mass flow, and vacuum thrust were 
calculated as described in the preceding section. The throat 
diameter was measured each test day to ensure that no distortion 
or eroding was occurring. None was observed, and an average 
value was used to calculate the throat area (tables I and II; one 
value for each throat section used). 

Efficiency Calculations . — The performance parameters 
(/ y, C h yy C*) were divided by the theoretical, one-dimen- 
sional-equilibrium (ODE) values obtained trom the Chemical 
Equilibrium Composition (CEC) program (rel. 21) to derive 
the efficiencies. The inlet enthalpy conditions were derived 
from measurements of the in jection pressure and temperature 
of the hydrogen and oxygen. Equations for the various ell icien- 
cies follow. The characteristic exhaust velocity efficiency is 


C* 


= — 

C *7 7; ( ODE ) 

the vacuum thrust coefficient efficiency 


is 


(9) 


C* 


F c,e^tSc 


m 


( 6 ) 






( 10 ) 


7. V.mODE) 


c - Fv 

l F.V “ - 


Pr,A 


f vHc 


l sp.v 


mg 


(7) 


( 8 ) 


and the vacuum specific impulse efficiency is 


sp, V 


' 'sp.v.n 


(ID 


sp. V'. Tin ODE ) 


O 

9 ; 

<D 

d 

a. 

o 

co 


1.02 i — 


0.98 


3 

CO 

V) 

E 

a 

"(5 

o 


o 

£ 0.94 
co 

-C 

o 

c 

o 

-3 

CO 

3 

n 

E 

Q 0.90 


o o 


Experimental data 

P c ,e/Pc,T = 0-994 - 0.00099 (O/F) 


l 


H V v 

Propellant mixture ratio, O/F 
Figure 1 3. — Effective chamber pressure correlation. 
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Heat Transfer 

Seven test firings were for experimental determinations of 
the attainable thrust performance of high-area-ratio rocket 
nozzles. During these firings, the outer wall temperatures of the 
carbon-steel, heat-sink nozzle skirt were measured. From these 
measurements, the heat fluxes were calculated. 

The technique of calculating heat flux from the wall tempera- 
ture data of a heat-sink nozzle was first employed in refer- 
ence 5. This technique, which is very simple in principle, 
resulted from observations that the nozzle wall temperature 
time histories were linear once the rocket thruster achieved 
steady-state thrust. The slope of the time response of the 
temperature was directly proportional to the heat flux to the 
wall. A detailed derivation of this principle is presented in 


appendix A of reference 5. Further analysis showed that axial 
conduction and radiation losses were the primary heat losses 
that would cause errors in the heat flux determination. In 
general, these were less than 2 percent in the area considered. 
This error analysis is detailed in appendix B of reference 5. 

During performance testing of the 1025: 1 nozzle, thermo- 
couple measurements were taken at a rate of 50 Hz, averaged 
in groups of five, and displayed at 0. 1 -sec intervals. Table IV 
lists the temperatures taken just prior to thruster shutdown. 
At this point, the thruster was at steady state with regard to 
the static pressure measurements in the nozzle. Measurements 
were taken at nine axial locations, which are identified in the 
table by the area ratio at the location. Also listed are the 
combustion conditions of the thruster for each of the firings. 


TABLE IV. — NOZZLE OUTSIDE WALL TEMPERATURES 


Reading 

Effective combustion 
chamber total pressure 
at nozzle entrance. 
p c.e 

Propellant 

mixture 

ratio, 

O/F 

MPa 

psi 

569 

12.326 

1787.7 

3.89 

570 

12.645 

1834.0 

5.97 

571 

12.488 

1811.1 

4.70 

575 

14.350 

2081.2 

4.65 

576 

14.605 

2118.2 

5.68 

577 

14.225 

2063.1 

4.47 

580 

16.364 

2373.3 

4.27 

601 

12.768 

1851.8 

6.15 

602 

12.542 

1819.0 

5.11 

603 

12.457 

1806.7 

4.01 


Expansion area ratio. € 
50 | 50.6 | 

Nozzle wall temperature 


K 

°R 

K 

°R 

K 

361.69 

651.05 

— 


325.56 

428.94 

772.09 


— 

370.09 

411.62 

740.91 

— 


363.23 

414.72 

746.50 



356.46 

438.81 

789.86 


— 

379.98 

416.22 

749.20 

— 

— 

375.08 

429.66 

773.38 

4)6.91 

750.44 

366.32 




421.42 

758.56 






406.07 

730.93 



°R 

586.00 

666.17 

653.82 

641.63 

683.96 

675.15 

659.37 


Reading 




E. 

<pansion area ratio 

£ 



^Z1 


101.2 

200 

202.4 

300 

303.6 | 





Nozzle wall 

temperature 





K 

°R 

K 

°R 

K 

°R 

K 

°R 

K 

°R 

569 



306.32 

55 1 .38 



299.35 

538.83 









570 

339.28 

610.71 

327.02 

588.63 







571 

337.71 

607.88 

327.87 

590.16 







575 

324.02 

583.23 

311.82 

561.27 







576 

577 

348.91 

353.87 

628.03 

636.96 

335.87 

345.15 

604.57 

621.27 













580 

330.98 

595.76 

317.27 

571.08 

365.07 

657.13 

328.47 

581.25 

317.13 

570.84 

601 









602 

391.28 

704.31 

365.24 

657.44 

350.88 

631.58 





603 

382.32 

688.18 

367.27 

661.08 

358.88 

645.99 






Reading 





Expansion 

area ratio. £ 




: 


388 

392.7 

500 

635 

800 

^ 1 


Nozzle wall temperature j 


K 

°R 

K 

°R 

K 

°R 

K 

°R 

K 

°R 

K ! 

°R 

569 

300.6 1 

541.09 


— 

297.56 

535.61 

301.44 

542.60 

303.82 

546.88 

302.86 

545.14 

570 

322.01 

579.61 



316.16 

569.09 

314.66 

566.38 

313.57 

564.42 

309.35 

556.83 

571 

324.53 

584.16 



319.16 

574.49 

319.45 

575.01 

320.53 

576.95 

317.% 

572.32 

575 

309.96 

557.92 

— 


303.44 

546.20 

304.79 

548.62 

304.29 

547.73 

304.32 

547.77 

576 

331.11 

595.99 



323.84 

582.92 

320.29 

576.52 

318.52 

573.34 

314.78 

566.61 

577 

Vp 49 

616 4Q 




600 95 

3V 06 

597 70 

330 83 

sqs so 

V>8 

590 90 

580 

313.56 

564.40 



307.35 

553.23 

308.32 

554.98 

306.74 

552.13 

307.64 

553.75 

601 

3 1 6.88 

570.38 

602 



347. 1 1 

624.80 









603 



358.74 

645.74 
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TABLED V — EXPERIMENTALLY MEASURED RATES OF INCREASE IN NOZZLE TEMPKRATURB 


Reading 

Effective combustion 
chamber total pressure 
at nozzle entrance, 

P,.e 

Propellant 

mixture 

ratio, 

O/F 

Expansion area ratio. £ 

50 

100 

200 | 

Measured rate of increase in nozzle temperature 

MPa 

psia 

K/s 

°R/s 

K/s 

°R/s 

K/s 

°R/s 

569 

12.326 

1787.7 

3.89 

63.67 

114.60 

34.89 

62.81 

18.76 

33.76 

570 

12.645 

1834.0 

5.97 

81.61 

146.89 

45.51 

81.91 

24.06 

43.31 

571 

12.488 

1811.1 

4.70 

73.04 

131.48 

40.52 

72.94 

21.29 

38.32 

575 

14.350 

2081.2 

4.65 

80.72 

145.29 

45.32 

81.57 

24.77 

44.59 

576 

14.605 

2118.2 

5.68 

88.77 

159.79 

49.77 

89.59 

26.85 

48.33 

577 

14.225 

2063.1 

4.47 

72.16 

129.88 

42.36 

76.25 

22.25 

40.05 

580 

16.364 

2373.3 

4.27 

82.68 

148.82 

48.15 

86.67 

26.69 

48.04 


Reading 

Expansion area ratio, £ 


300 

.488 


500 

635 

800 

975 


Measured rate of increase in nozzle temperature I 


K/s 

°R/s 

K/s 

°R/s 

K/s 

°R/s 

K/s 

°R/s 

K/s 

°R/s 

K/s 

°R/s 

569 

13.56 

24.40 

10.92 

19.65 

8.26 

14.86 

7.00 

12.60 

6.57 

11.82 

5.60 

10.08 

570 

16.89 

30.40 

13.37 

24.06 

9.93 

17.87 

7.69 

13.84 

6.03 

10.86 

5.03 

9.06 

571 

15.24 

27.44 

12.13 

21.84 

9.13 

16.44 

7.09 

12.77 

5.76 

10.37 

4.73 

8.52 

575 

17.45 

31.41 

13.74 

24.74 

10.26 

18.46 

8.12 

14.62 

6.41 

1 1 .53 

5.43 

9.78 

576 

18.63 

33.54 

14.77 

26.58 

10.99 

19.79 

8.58 

15.45 

6.59 

11.87 

5.63 

10.13 

577 

16.03 

28.85 

12.73 

22.92 

9.53 

17.15 

7.58 

13.64 

6.24 

11.23 

5.41 

9.73 

580 

18.71 

33.68 

14.91 

26.83 

1 1.21 

20.17 

8.64 

15.55 

6.79 

12.22 

5.83 

10.49 


The time rate of change of the nozzle wall temperature 
measurements was also noted. For every thermocouple, the rate 
of change for any thermocouple was constant during the last 
second of firing. These measurements are tabulated in table V, 
which shows the rate of temperature increases dT/dt , for all 
nine locations for each of seven firings. From these values of 
dT/dt. we could calculate the temperature of the inside wall: 


l 2 l 2 

Q= | c /" ()A s = J 
L i /- 1 


c£nl)i 
cos 0 , 


(14) 


Calculated total heat rates for the seven firings are tabulated in 
table VII. 


ML 

4a dt 


1- 




\Koj 


In 


v j 


( 12 ) 


Equation (12) is derived in reference 5. Values for T () were 
obtained from table IV, and values for dT/dt were obtained 
from table V. Calculated inside wall temperatures are pre- 
sented in table III. The heat fluxes to the wall of the nozzle were 
also calculated: 




ft 


kRi dT 
2 a dt 


{ 


V 


K 

*i 


-1 


(13) 


Equation (13) is also derived in reference 5. Calculated nozzle 
wall heat fluxes are tabulated in table VI. 

The heat rate Q to the walls of a rocket nozzle between two 
axial locations can be determined by integrating the heat flux 
values with respect to the nozzle surface area. Details of this 
derivation are in reference 5: 


Boundary Layer 

Two different instrumentation configurations were used to 
obtain total pressure measurements in the 12 firings reported 
herein. The first was with a single, nominally 4-in. high, 
boundary-layer rake (fig. 10(a)) with seven total pressure 
probes. The other was with two small rakes (figs. 1 0(b) and (c )), 
nominally 1-in. high with three and four total pressure probes, 
respectively. The rakes were mounted at the exit of the nozzle, 
with the two short rakes mounted 1 80° apart. This arrangement 
allowed total pressure measurements to be obtained at 14 dif- 
ferent distances from the nozzle wall. The rakes were installed 
with the probe tubes inserted into the exit of the nozzle, with the 
probe tubes parallel to the nozzle wall. Because of the nozzle 
wall divergence, the body of the rake was not radial, but per- 
pendicular, to the wall. As a result, the individual probes were 
each located at a slightly different axial dimension, and hence, 
at a slightly different expansion area ratio (see f igs. 1 0(a) to (c )). 

The boundary layer rake data were not manipulated, except 
for the total pressure measurements, which were normalized by 
the effective chamber pressure. 
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TABLE VI— EXPERIMENTALLY DETERMINED NOZZLE HEAT FLUX DISTRIBUTION 


Reading 

Effective combustion 
chamber total pressure 
at nozzle entrance, 

Pc.e 

Propellant 
mixture ratio, 
O/F 

Characteristic exhaust 
velocity efficiency, 
(best fit curve) 

percent 

Expansion area ratio, £ 

50 

100 

Heat flux to nozzle walls, as measured 

MPa 

psia 

kW/m : 

Btu/in.“-sec 



569 

12.326 

1787.7 

3.89 

99.82 

1422.77 

0.8706 

772.02 

0.4724 

570 

12.645 

1834.0 

5.97 

99.07 

1823.66 

1.1 159 

1006.70 

.6160 

571 

12.488 

1811.1 

4.70 

99.65 

1632.29 

.9988 

896.55 

.5486 

575 

14.350 

2081.2 

4.65 

99.67 

1803.88 

1.1038 

1002.61 

.6135 

576 

14.605 

2118.2 

5.68 

99.26 

1983.81 

1.2139 

1101.16 

.6738 

577 

14.225 

2063.1 

4.47 

99.72 

1612.51 

.9867 

937.24 

.5735 

5 SO 

16.364 

2373.3 

4.27 

99.76 

1847.68 

1.1306 

1065.20 

.65 1 8 


Reading 

Expansion area ratio, e | 



200 


300 

388 

5™ | 


Heat flux to nozzle walls, as measured 


kW/m : 

Btu/in.~-sec 

kW/m 2 

Btu/in.“-sec 

kW/irf 

Btu/in.'-sec 



569 

41 1.99 

0.2521 

296.78 

0.1816 

238.60 

0.1460 

180.26 

0.1103 

570 

528.52 

.3234 

371.30 

2272 

292.20 

.1788 

216.70 

.1326 

571 

467.56 

.2861 

333.71 

.2042 

265.24 

.1623 

199.38 

.1220 

575 

544.20 

.3330 

382.09 

.2338 

300.37 

.1838 

223.89 

.1370 

576 

589.80 

.3609 

408.07 

.2497 

322.76 

.1975 

239.91 

.1468 

577 

488.80 

.2991 

350.87 

.2147 

278.31 

.1703 

207.88 

.1272 

580 

586.20 

.3587 

409.71 

.2507 

325.87 

.1994 

244.65 

.1497 


Reading 

Expansion area ratio, e j 


635 

800 


975 i 


Heat 11 ux to nozzle walls, as measured 





7 




kW/nf 

Btu/in.“-sec 

kW/nT 

Btu/in.'-sec 

kW/m' 

Btu/in.”-see 

569 

152.64 

0.0934 

143.00 

0.0875 

121.91 

0.0746 

570 

167.67 

.1026 

131.39 

.0804 

109.49 

.0670 

571 

154.60 

.0946 

125.51 

.0768 

102.96 

.0630 

575 

176.99 

.1083 

139.40 

.0853 

118.16 

.0723 

576 

187.12 

.1145 

143.65 

.0879 

122.41 

.0749 

577 

165.22 

.101 1 

135.81 

.0831 

117.67 

.0720 

580 

188.27 

.1152 

147.90 

.0905 

126.82 

.0776 
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TABLE VII.— TOTAL HKAT RATH VALUES ADJUSTED TO A COMMON P, g 


Reading 

Effective 

combustion 

Propellant 
mixture ratio. 

Heat rate from £ = 1 40 to 1 025 


chamber total 

O/F 

Experimental 


pressure at 
nozzle entrance. 

P .... 


No adjustments 

Adjusted to 
P ce = 2063.1 psia 

Ad justed to 
r|c< =100 percent 

Adjusted to 
* = 100 percent and 











P ce = 2063.1 psia 


MPa 

psia 


kW 

Btu/s 

kW 

Btu/s 

kW 

Btu/s 

kW 

Btu/s 

569 

12.326 

1787.7 

3.89 

380.43 

360.67 

426.63 

404.47 

381.81 

361.97 

428.17 

405.93 

570 

12.645 

1834.0 

5.97 

425.68 

403.56 

467.71 

443.41 

433.71 

411.18 

476.54 

451.78 

571 

12.488 

1811. 1 

4.70 

390.23 

369.96 

433.10 

410.60 

w 

NO 

TO 

'o 

oc 

372.56 

436.15 

413.49 

575 

14.350 

2081.2 

4.65 

443.65 

420.60 

440.56 

417.67 

446.59 

423.39 

443.48 

420.44 

576 

14.605 

21 18.2 

5.68 

471.11 

446.63 

461.27 

437.31 

478.15 

453.31 

468.17 

443.85 

577 

14.225 

2063.1 

4.47 

416.83 

395.17 

416.83 

395.17 

419.17 

397.39 

419.17 

397.39 

580 

16.364 

2373.3 

4.27 

476.84 

452.07 

426.30 

404.15 

479.14 

454.25 

428.35 

406.10 


Reading 

Heat rate from £ = 140 to 1025 

TDK/MABL turbulent 

Adjusted to 
r|< 7 * =100 percent 

Adjusted to 
rif'. = 100 percent and 
P c e = 2063. 1 psia 

kW 

Btu/s 

kW 

Btu/s 

569 

520.36 

493.33 

583.57 

553.25 

570 

617.50 

585.42 

678.48 

643.23 

571 

559.67 

530.59 

621.14 

588.87 

575 

626.02 

593.50 

621.66 

589.36 

576 

681.88 

646.45 

667.65 

632.96 

577 

608.06 

576.47 

608.06 

576.47 

580 

675.04 

639.97 

603.48 

572.13 


Results 

Performance Results 

Atmospheric Pressure Tests . — Tests were performed at 
atmospheric pressure to determine the relationship between the 
ef fective and measured chamber pressures of the thruster. The 
tests were conducted with low-arca-ratio configurations 
(e = 10.7:1 and 4:1), the performance of which is well doc- 
umented and agrees with calculated values from the TDK 
program. Because of the low-area-ratio of the nozzles, an 
altitude condition was not necessary to obtain full, unseparated 
nozzle flow. The results of the nine successful atmospheric 
tests are summarized in table VIII. In this table, the measured 
combustion chamber static pressure at the injector lace is listed 
as P , and equation (4) was used to derive the P ( T values from 
the P values. The effective chamber pressures P ( . € , derived 
from thrust measurements as previously described, are also 
listed in table VIII. A consistent variation between P c e and P c J 
was observed and was attributed to variations in the static 
pressure profile that most likely occurred at the static tap used 
for the P ca measurements. 

So that the decrease in thrust attributable to combustion 
losses could be properly accounted for, the characteristic exhaust 
velocity C* and the characteristic exhaust velocity efficiency 


r\^, were derived for both the atmospheric and altitude tests. 
Within the range of these tests, chamber pressure had no el feet 
on r\ c * and only a slight variation with respect to OIF. Fig- 
ure 14 shows r| c * as a function of OIF for all atmospheric and 
altitude firings. A mean value of r \ c * was described by a 
second-order polynomial curve fit (eq. (15)) by the least- 
squares method, with values ranging from approximately 
99.0 to 99.9 percent. 

lie* =98.43 + 0.824 (O/F)- 0.1 20 (O/F) 2 05) 

Altitude Tests. — High-area-ratio nozzle tests were performed 
at altitude conditions to avoid separated flow in the divergent 
portion of the nozzle. The first test objective was to ascertain 
whether the flow was attached or separated by examining the 
nozzle wall static pressure distribution. Static pressures, which 
were measured at eight axial locations, are given in table IX. 
Figure 15 shows a typical distribution along the length of the 
nozzle. Plotted there from reading 577 are the static pressure 
ratios P s /P. e versus the nozzle expansion ratio of the pressure 
tap locations. When plotted on log-log coordinates, the result is 
a straight line. If the flow were separated, the pressure distribu- 
tion would display a sudden increase. As this was not the case 
for any of the tests, all the data reported have attached flow . 
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TABLE VIII — RESULTS OF ATMOSPHERIC PRESSURE TESTS 


Reading 

Expansion 
area ratio, 
e 

Nozzle throat 
area, 

A t 

Measured chamber pressure 

Propellant 

mixture 

ratio. 

O/F 

Measured 

thrust, 

F 

At injector face, 
P c.A 

Corrected for momentum 
pressure loss, 

P c,T 

cni~ 

in 

2 

MPa 

psia 

MPa 

psia 

N 

lb f 

514 

10.72 

5.103 

0.7909 

13.942 

2022.1 

13.893 

2014.9 

4.21 

11 209 

2520.0 

515 

10.72 

5.103 

.7909 

15.801 

2291.7 

15.748 

2284.0 

4.16 

12 677 

2850. 1 

523 

3.99 

5.091 

.7890 

12.254 

1777.3 

12.211 

1771.0 

4.19 

9491 

2133.7 

524 

3.99 

5.091 

.7890 

12.524 

1816.4 

12.461 

1807.3 

5.38 

9 735 

2188.6 

526 

4.02 

5.047 

.7823 

14.362 

2083.0 

14.293 

2072.9 

5.66 

11 174 

2512.2 

527 







14.746 

2138.7 

14.675 

2128.4 

5.60 

I 1 285 

2537.1 

528 







15.096 

2189.4 

15.023 

2178.8 

5.50 

11511 

2587.8 

529 







12.825 

1860.1 

12.756 

1850.1 

5.83 

9 808 

2205.0 

530 







14.642 

2123.6 

14.564 

2112.3 

5.98 

1 1 423 

2568.2 


Reading 

Vacuum thrust, a 

p v 

Propellant flow 
rate, 
rh 

Fuel injection 
pressure. 

p « 

Fuel injection 
temperature. 

T fi 

Oxidizer 

injection 

pressure, 

P oi 1 


N 

'bf 

kg/s 

,b n/ s 

MPa 

psia 

K 

°R 

MPa 

psia 

514 

1 1 746 

2640.7 

2.852 

6.287 

17.818 

2584.2 

294.9 

530.8 

15.271 

2214.8 

515 

13 214 

2970.8 

3.158 

6.962 

20.030 

2905.0 

294.7 

530.4 

17.411 

2525.1 

523 

9 690 

2178.5 

2.5(X) 

5.512 

16.147 

2341.9 

308.5 

555.3 

13.399 

1943.3 

524 

9 934 

2233.3 

2.644 

5.828 

15.389 

2231.9 

306.5 

551.7 

13.891 

2014.6 

526 

1 l 374 

2557.0 

3.037 

6.6% 

17.508 

2539.3 

309.3 

556.7 

16.326 

2367.8 

527 

1 1 486 

2582.2 

3.102 

6.839 

17.915 

2598.3 

300.8 

541.5 

16.791 

2435.3 

528 

11711 

2632.8 

3.143 

6.928 

18.258 

2648.0 

299.1 

538.4 

17.036 

2470.8 

529 

10 008 

2250.0 

2.707 

5.968 

15.341 

2224.9 

299.8 

539.6 

14.329 

2078.2 

530 

11 424 

2568.3 

3.082 

6.794 

17.420 

2526.4 

300.9 

541.7 

16.602 

2407.8 


Reading 

Oxidizer 

injection 

temperature, 

T 

oi 

Theoretically predicted 

Effective chamber 
total pressure 
calculated from 
thrust. 

p c.e 

Correlation 
pressure ratio 
for use in 
altitude tests. 

p c.c /P c.T 

Characteristic 

exhaust 

velocity 

efficiency, 

fic*- 

percent 

ODE vacuum 
thrust 

coefficient. 

C F.V,Th 

(ODE) 

TDK vacuum 
thrust 

coefficient, 

C F.V.Th 

(TDK) 

Vacuum 

thrust 

coefficient 

efficiency. 

TDK/ODE, 

V 

Cp,V 

K 

°R 

MPa 

psia 

514 

90.9 

163.6 

1.737 

1.673 

96.3 

13.758 

1995,4 

0.990 

99.1 

515 

88.8 

159.9 

1.736 

1.673 

%.4 

15.483 

2245.5 

.983 

100.6 

523 

%.7 

174.1 

1.601 

1.557 

97.2 

12.230 

1773.8 

!.(X)2 

100.2 

524 

92.9 

167.3 

1.615 

1.570 

97.2 

12.433 

1803.2 

.998 

99.5 

526 

110.2 

198.3 

1.619 

1.573 

97.2 

14.327 

2077.9 

1 .002 

99.8 

527 

110.2 

198.4 

1.619 

1.572 

97.1 

14.476 

2099.5 

.986 

98.5 

528 

92.3 

166.2 

1.617 

1.571 

97.2 

14.769 

2142.0 

.983 

98.9 

529 

92.8 

167.0 

1.622 

1.576 

97.2 

12.585 

1825.3 

.987 

99.0 

530 

93.1 

167.6 

1.623 

1.576 

97.1 

14.365 

2083.3 

.986 

99.6 


Reading 

Vacuum 

specific 

impulse 

[ sp.V' 

s 

Vacuum 

specific 

impulse 

efficiency. 

Xp.V’ 

percent 

Ambient pressure 
around nozzle. 

P a 

kPa 

psia 

514 

420.0 

95.5 

98.143 

14.234 

515 

426.7 

%.9 

98.109 

14.229 

523 

395.2 

97.4 

| 

i 

524 

383.2 

%.7 

1 

I 

526 

381.9 

%.9 

♦ 

♦ 

527 

377.6 

95.7 

98.854 

14.337 

528 

380.0 

96.0 

98.819 

14.332 

529 

377.0 

%.2 

98.785 

14.327 

530 

378.0 

96.8 

98.681 

14.312 


Pleasured thrust corrected to vacuum conditions. 
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Static pressure ratio, P s /P c>e 31 Characteristic exhaust velocity efficiency, t\ c 



3 4 5 6 


Propellant mixture ratio, O/F 

gure 14. — Characteristic exhaust velocity efficiency as function of propellant mixture 
ratio. 



Expansion area ratio, e 


Figure 15. — Typical nozzle wall static pressure distribution (reading 577). 
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TABLE IX — NOZZLE WALL STATIC PRESSURES 


Reading 

Effective combustion at 
chamber total pressure 
nozzle entrance. 
P c.e 

Propellent 

mixture 

ratio, 

O/E 


Expansion 

area ratio, e 

1 

100 

101.2 I 

200 | 

Nozzle wall static pressure. P s 


MPa 

psia 


Mvjm 


Wfi ?■ 

KH 



569 

12.326 

1787.7 

3.89 

13.34 

1.935 


— 

5.766 

0.8362 

570 

12.645 

1834.0 

5.97 

14.39 

2.087 

— 

— 

6.281 

.9109 

571 

12.488 

1811.1 

4.70 

14.90 

2.161 

— 

— 

6.470 

.9383 

575 

14.350 

2081.2 

4.65 

17.03 

2.470 

— 

— 

7.350 

1.066 

576 

14.605 

2118.2 

5.68 

17.20 

2.495 

— 

— 

7.426 

1.077 

577 

14.225 

2063.1 

4.47 

16.80 

2.436 

— 

— 

7.302 

1.059 

580 

16.364 

2373.3 

4.27 

19.09 

2.769 

— 

— 

8.253 

1.197 

601 

12.768 

1851.8 

6.15 

— 

— 

14.78 

2.143 

— 


602 

12.542 

1819.0 

5.11 

— 

— 

14.79 

2.145 

— 

— 

603 

12.457 

1806.7 

4.01 

— 

— 

13.68 

1.984 

— 

— 


Reading 




Ex 

pansion 

area ratio, e 





202.4 

3tX) 

303.6 

388 

392.7 | 


Nozzle wall static pressure, P s 







bseh 


B350BI 

kPa 

MBEM 

569 


— 

3.476 

0.504 1 


— 

2 S22 

0.3658 

— 

— 

570 


— 

3.929 

.5699 

— 

— 

2.990 

.4337 

— 

— 

571 

■ 


3.895 

.5649 


— 

2.832 

.4108 

— 


575 

E 


4.456 

.6462 

— 


159 

.4717 

— 


576 

— 


4.656 

.6753 


— 

119 

.5093 

— 

— 

577 

— 

— 

4.410 

.6396 

— 

— 

Esa 

.4688 

— 


580 

— 

— 

4.955 

.7186 

— 

— 

3.609 

.5234 

— 

— 


6.847 

0.9930 

— 

— 

4.028 

0.5842 

— 

— 

3.026 

0.4389 


6.723 

.9750 

— 

— 

4.003 

.5805 

— 

— 

2.968 

.4305 


6.172 

.8952 

— 

— 

3.725 

.5403 

— 

— 

2.755 

.3996 


Reading 



Expansion 

area ratio, e 


1 

■ 

500 

635 

800 

975 | 


Nozzle wall static pressure, P s 

■ 

■3ZI 

KSEB 

non 

■389 

tm 


kPa 

■SEE 

569 

1.789 

0.2594 

1.351 

0.1959 

1 .008 

0.1462 

0.7853 

0.1139 

570 

2.224 

.3225 

1.624 

.2356 

1,197 

.1736 

.9163 

.1329 

571 

2.035 

.2952 

1 .496 

.2169 

1.105 

.1602 

.8550 

.1240 

575 

2.299 

.3335 

1.687 

.2446 

1.247 

.1809 

.9646 

.1399 

576 

2.535 

.3676 

1.854 

.2689 

1.362 

mm 

um 

.1520 

577 

2.289 

.3320 

1.674 

.2428 

1 .246 

mm 

.9550 

.1385 

580 

2.549 

.3697 

1.872 

.2715 

1.380 

■i 

1.069 

.1551 




— 

— 


be 

— 

— 


— 


— 



m 

— 

— 


— 


— 

— 


■ 

— 

— 


Ten successful firings were accomplished at altitude — seven 
with the 1025: 1 -area-ratio nozzle and three with the nozzle 
truncated to an area ratio of 440:1. Table II summarizes the 
hoi- lire results, including measured and calculated values. 

Figure 16 shows the nozzle thrust performance in terms of 
C j. v . Two sets of data are shown: the first is for the original 
nozzle with the 1025:1 area ratio, and the second is for the 
truncated nozzle with the 440: 1 area ratio. Straight lines of the 


best fit by the least-squares method are shown. For the 1025: l 
nozzle, the thrust coefficients ranged from approximately 
1 .92 to 2.02. and for the 440: 1 nozzle, they ranged from 1.83 
to 1.94. 

The nozzle thrust efficiency is shown in figure 1 7 as straight 
lines of the best fit by the least-squares method. The efficiencies 
ranged from approximately 96.6 to 97.5 percent for the 1 025: 1 
nozzle and from 94.0 to 94.2 percent for the 440:1 nozzle. 
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Propellant mixture ratio, O/F 


Figure 16. — Nozzle thrust performance as a function of mixture ratio. 



Propellant mixture ratio, O/F 

Figure 1 7.— Nozzle thrust efficiency as a function of mixture ratio. 
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Figure 18. — Overall thruster performance as a function of mixture ratio. 
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Figure 1 9. — Overall thruster efficiency as a function of mixture ratio. 


Figure 1 8 shows the overall thruster performance with a plot 
of specific impulse versus O/F for both the 1025: 1 and 440: 1 
configurations. The faired curves were obtained from the 
product of the faired curves of figure 17 and the theoreti- 
cal ODE values of reference 2 1 . The specific impulse attained 
was as high as 488 sec for the 1 025: 1 -nozzle configuration and 
467 sec for the 440: 1 -nozzle configuration. 


Figure 19 shows the overall thruster efficiency as the spe- 
cific impulse efficiency plotted as a function of O/F for the 
1025:1- and 440: 1 -area-ratio configurations. Ideally, overall 
specific impulse efficiency should be equal to the product of 
r \ c * and rj^ ^as shown in equation ( 16): 

X Tl/ ViV (16) 
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The faired curves shown in figure 19 were obtained from 
the product of the best-fit curves of T| c * and T| c v/ from figures 
14 and 17, respectively. The coincidence of the faired curves 
through the center of the apparent data scatter reinforces the 
quality of the results. Deviations from this relationship were 


attributed to measurement uncertainties in effective combus- 
tion chamber total pressure, vacuum force, and mass flow rate. 
Values of r| ; ^ ^ranged from 95.5 to 97.5 percent for the 1025:1 
configuration and from 93.3 to 94.0 percent for the 440:1 
configuration. 


TABLE X. — TDK/BLM TURBULENT PREDICTIONS 


Reading 

Nozzle exit 
expansion 
area ratio, 
e 

Effective combustion 
chamber total pressure 
at nozzle entrance. 

I’m 

Measured 
propellant 
mixture ratio, 
O/F 

Predicted propellant 
flow rate 

MPa 

psia 

k.g/s 

tt>m/s 

569 

1025 

12.326 

1787.7 

3.89 

2.5144 

5.5432 

570 



12.645 

1 834.0 

5.97 

2.7306 

6.0198 

571 



12.488 

1811.1 

4.70 

2.5937 

5.7181 

575 



14.350 

2081.2 

4.65 

2.9469 

6.5628 

576 



14.605 

21 18.2 

5.68 

3.1214 

6.8813 

577 



14.225 

2063. 1 

4.47 

2.9370 

6.4748 

580 



16.364 

2373.3 

4.27 

3.3617 

7.4 1 1 1 

601 

440 

12.768 

1851.8 

6.15 

2.7384 

6.0371 

602 

440 

12.542 

1819.0 

5.1 1 

2.6018 

5.7359 

603 

440 

12.457 

1 806.7 

4.01 

2.5135 

5.5413 


Reading 

Nozzle exit 
expansion 
area ratio, 
e 

Computer code 

TDK/BLM Turbulent 

Predicted characteristic 
exhaust velocity, 

Predicted vacuum thrust. 

P v 

Predicted 

vacuum 

thrust 

coefficient. 

C F.V 

Predicted vacuum 
thin st coefficient 
efficiency, 

Vv- 

percent 

m/s 

ft/s 

N 

ih, 

569 

1025 

2483.97 

8149.51 

l 1 701 

2630.72 

1.87 

95.49 

570 



2346.55 

7698.65 

12 754 

2867.32 

1.99 

94.84 

571 



2439.50 

8003.61 

12 132 

2727.59 

1.92 

95.26 

575 



2442.52 

8013.52 

1 3 935 

3132.98 

1.92 

95.45 

576 



2370.80 

7778.23 

14619 

3286.56 

1.98 

95.12 

577 



2454.20 

8051.85 

1 3 742 

3089.58 

1.91 

95.48 

580 



2466.44 

8091.99 

15 719 

3533.93 

1.90 

95.56 

601 

440 

2334.25 

7658.31 

12 445 

2797.95 

1.95 

93.61 

602 

440 

2413.37 

7917.89 

1 1 861 

2666.66 

1.89 

93.75 

603 

440 

2481.10 

8140.08 

1 1 410 

2565.17 

1.83 

93.83 


Reading 

Nozzle exit 
expansion 
area ratio. 
£ 

Computer code 

Predicted 

vacuum 

specific 

impulse 

(adjusted), 

hp.V 

Predicted vacuum 
specific impulse 
efficiency (adjusted), 

T 

UpV’ 

percent 

ODE 

ODK 

MOC 

TDK/BLM 

turbulent 

Predicted vacuum specific impulse. 

Up.V’ s 

569 

1025 

500.63 

499.61 

495. 1 3 

474.59 

473.73 

94.63 

570 



505.53 

502.12 

498.63 

476.31 

471.89 

93.35 

571 



504.43 

502.71 

498.57 

477.01 

475.35 

94.23 

575 



503.90 

502.51 

498.30 

477.38 

475.79 

94.42 

576 



505.55 

503.10 

499.40 

477.61 

473.96 

93.75 

577 



503.48 

502.22 

497.93 

477.17 

475.82 

94.51 

580 



502.54 

501.62 

497.23 

476.84 

475.70 

94.66 

601 

440 

497.81 

494.39 

482.30 

463.46 

458.64 

92.13 

602 

440 

498.99 

496.98 

483.92 

464.91 

462.62 

92.71 

603 

440 

496.58 

495.52 

481.57 

462.92 

462.01 

93.04 
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TABLE XI — TDK/M ABL TURBULENT PREDICTIONS 


Reading 

Nozzle exit 
expansion area 
ratio, 
e 

Effective 
combustion 
chamber total 
pressure at nozzle 
entrance. 

Pcx- 

Propellant 

mixture 

ratio. 

O/F 

Predicted propellant 
flow rate 

MPa 

psia 

kg/s 


569 

1025 

12.326 

1787.7 

3.89 

2.5095 


570 



12.645 

1834.0 

5.97 

2.7240 


571 



12.488 

1811.1 

4.70 

2.5883 


575 



14.350 

2081.2 

4.65 

2.9708 


576 



14.605 

21 18.2 

5.68 

3.1143 

6.8658 

577 



14.225 

2063.1 

4.47 

2.9311 

6.4618 

580 



16.364 

2373.3 

4.27 

3.3556 

7.3976 

601 

440 

12.768 

1851.8 

6.15 

2.7342 

6.0277 

mm 

440 

12.542 

1819.0 

5.1 1 

2.5989 

5.7296 

m 

440 

12.457 

1806.7 

4.01 

2.5127 

5.5395 


Reading 

Nozzle exit 
expansion area 
ratio, 

£ 

Predicted 
characteristic 
exhaust velocity. 
C* 

Computer code 

TDK/MABL Turbulent 

Predicted vacuum 
thrust, 

Pv 

Predicted 

vacuum 

thrust 

coefficient, 

C F.V 

Predicted vacuum 
thrust coefficient 
efficiency, 

v ■ 

MV 

percent 

m/s 

ft/s 

N 

ib, 

569 


2488.86 

8165.54 

masM 

2628.16 

1.87 

95.40 

570 


■ 

2352.27 

7717.43 

mam 

2866.68 

1.99 

94.82 

571 



2444.64 

8020.47 

UB 

2^22 2 1 

1.91 

95.07 

575 



2447.53 

8029.95 

B 1 

3130.1 1 

1.91 

95.36 

576 



2376.16 

7795.79 

BIS 

3285.59 

1.98 

95.09 

577 



2459.12 

8067.99 

BB 

3087.30 

1.91 

95.41 

580 



2470.97 

8106.85 

BH 

3531.83 

1.89 

95.50 

601 


2337.88 

7670.21 

BB 

2800.75 

1.95 

93.70 

602 


2416.01 

7926.55 

BB 

2671.58 

1.89 

93.93 

603 


2481.91 

8142.74 

UB9 

2568.73 

1.83 

93.96 


Reading 

Nozzle exit 
expansion area 
ratio, 
e 

Computer code | 

ODE 

ODK 

MOC 

TDK/ 

MABL 

turbulent 

Predicted 

vacuum 

specific 

impulse 

(adjusted), 

Up.v 

Predicted vacuum 
specific impulse 
efficiency 
(adjusted). 

’1, 

Lp.V’ 

perceni 

Predi 

cted vacuum specific impulse, 
Isp.V. s 

569 

1025 

500.63 

499.61 

495.13 

475.06 

474.20 

94.72 

570 



505.53 

502.12 

498.63 

477.37 

472.93 

93.55 

571 



504.43 

502.71 

498.58 

477.07 

475.41 

94.25 

575 



503.90 

502.51 

498.32 

477.92 

476.33 

94.53 

576 



505.55 

503.09 

499.43 

478.54 

474.89 

93.94 

577 



503.48 

502.22 

497.95 

477.78 

476.42 

94.63 

580 



502.54 

501.62 

497.24 

477.43 

476.29 

94.78 

601 

440 

497.81 

494.36 

482.36 

464.64 

459.82 

92.37 

602 

440 

498.99 

496.98 

483.98 

466.28 

463.98 

92.98 

603 

440 

496.58 

495.52 

481.61 

463.71 

462.80 

93.20 


All the results discussed previously were compared with 
analytical predictions obtained from the TDK computer pro- 
gram, as previously described. Two turbulent models were run 
for each firing: one with the BLM module and one with the 
MABL module. A laminar boundary layer module was also run 


with the MABL module. Tables X and XI give the turbulent 
results for the BLM and MABL modules.(Laminar results are 
given in table XVI in app. B.) Close examination of these 
tabulated results shows that there is no significant difference 
between the results obtained with the turbulent BLM and 
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MABL modules. This is consistent with the findings of refer- 
ence 6, which studied a 300:1 liquid hydrogen and liquid 
oxygen nozzle. For simplicity, the rest of this discussion is 
limited to the results obtained with the MABL module. 

Table II gives the individual firing data points with the scatter, 
and the following paragraph discusses the mean values repre- 
sented by curve fits of these data. Figure 20 is a plot of the pre- 
dicted results and the attained I yV of the thruster with the 1025: 1 
nozzle as a function of Off. Tne ODE values are the predicted 
ideal, one-dimensional equilibrium values of specific impulse. 
The ODK values are the predicted results for one-dimensional, 
nonequilibrium flow, and the drop in v from ODE to ODK 
represents the predicted loss in performance due to kinetics. 
The third line, which was obtained from the MOC module, 
represents the inviscid. two-dimensional, nonequilibrium pre- 
dictions. The difference between the ODK and MOC values is 
the loss in performance due to nozzle divergence and shock 
losses occurring in the inviscid core portion of the nozzle flow. 

The next specific impulse decrement to be considered is the 
losses attributable to the boundary layer. The analytical model 


used for this is the MABL module of the TDK computer 
program. The first line, which is labeled TDK/M ABL(lam), 
represents losses attributable to laminar boundary layer growth 
along the nozzle wall. This is followed by the line labeled 
“laminar/' which represents performance losses attributable 
to combustion losses or energy release losses. This was 
determined by multiplying the last predicted / v values by 
r| c *. These predicted / v , values can be compared with the 
experimentally attained / v . Next, is the line representing the 
experimental results. For simplicity, the figure — the individual 
scatter of the experimental results — was not included and only 
the best fit curve was shown. The line labeled TDK/MABL 
(turb) represents losses attributable to turbulent boundary layer 
growth along the nozzle wall. The last performance decrement 
to be considered is the / v losses attributable to combustion 
losses or energy release loss. These are given in the bottom line 
on the chart labeled “turbulent." The experimentally measured 
values did not fall on either the predicted laminar values or 
turbulent values as expected, but do have very similar shapes 
and are very nearly parallel to one another. 



Propellant mixture ratio, O/F 


Figure 20. — Predicted thrust chamber losses from ideal performance. Area ratio, e, 
1025:1. 
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Figure 21 is a plot of \] { as a function of O/F. Values values, and turbulent values are uniformly 2-percent lower than 

of the predicted laminar and turbulent impulse efficiencies are the experimentally achieved values. 

shown for comparison to the curve of the experimentally Figure 22 is a plot of Cy v for the 1025:1 nozzle as a function 

achieved impulse efficiency. Predicted laminar values are of OIF . Shown here are the values of the predicted laminar, 

about 0. 5-percent higher than the experimentally achieved turbulent, and experimentally achieved nozzle C F y. For 



3 4 5 6 7 


Propellant mixture ratio, O/F 

Figure 21 . — Thrust chamber performance efficiency. Area ratio, e, 1025:1 . Specific 
impulse, r\. , is based on ideal one-dimensional equilibrium (ODE) results. 

l sp,v 




Propellant mixture ratio, O/F 


Figure 22. — Experimental and predicted nozzle vacuum thrust coefficient. Area ratio, e 
1025:1. 
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Vacuum specific impulse, l sPiV , s Vacuum thrust coefficient efficiency, r,^ percent 



Propellant mixture ratio, OIF 


Figure 23. — Nozzle thrust coefficient efficiency. Area ratio, £, 1025:1 . 



Propellant mixture ratio, OIF 


Figure 24. — Thrust chamber performance efficiency. Area ratio, e, 440:1 . 
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clarity, the best fit straight line was plotted instead of the 
experimental data scatter. As in figures 20 and 21, the experi- 
mentally achieved results are greater than the values obtained 
from the turbulent analysis but only slightly less than the values 
obtained from the laminar analysis. 

Figure 23 shows nozzle X] Cfr ^ as a function of OIF. As in 
figures 2 1 to 22. the experimentally achieved values are about 
0.5-percent lower than the predicted laminar values and uni- 
formly 2-percent higher than the predicted turbulent values. 

For the next series of comparisons, the rearmost part of the 
nozzle was removed. This provided a nonoptimized nozzle 
truncated at an exit area ratio of 440: 1 instead of extending to 
the optimized 1025:1 area ratio. The resulting divergence 
angles were higher than they would have been for an optimized 
440: 1 configuration. Figure 24 is a plot of predicted thruster 
performance and attained thruster performance for the 440:1 
configuration. Shown are the results of three firings. In compar- 
ing the 440:1 performance (fig. 24) to that of the 1025:1 
configuration (fig. 20), one observes the obvious decrease in 
performance (~20-s / ) attributable to both a reduced area ratio 
and an increased nozzle exit divergence angle. In comparison, 
the attained performance represented by the best fit curve to the 
predicted performance for the 440: 1 configuration falls much 
closer to the predicted turbulent values than does the 1025:1 
configuration, and it is lower than the predicted laminar values. 
Figure 25 summarizes this same result by showing r|j ^ for the 
440:1 configuration. In comparison to the 1025:1 configura- 


tion, again the experimentally attained efficiency for the 440: 1 
configuration is much closer to the TDK predicted turbulent 
values than to the predicted laminar values. 

Figure 26 shows the C F y of the 440: 1 -nozzle configuration 
as a function of OIF. As was the case for the 1025: 1 configura- 
tion, C F y varies linearly with respect to OIF , and the experi- 
mental values fall quite close to the TDK turbulent prediction, 
and well below the linear predictions. Again, this is signifi- 
cantly different from the 1025:1 results, where the experi- 
mental values were closer to the laminar predictions. This 
correspondence is further illustrated in figure 27, which is a 
plot of V[ Cfr ^ as a function of OIF. The efficiency expressed 
here is, as elsewhere in this paper, based on the ODE values. 
The experimentally attained efficiency is about 0.25-pereeni 
higher than the turbulent prediction and nearly 2-percent lower 
than the laminar predictions. 

Heat Transfer Results 

The following is a discussion of the heat-transfer results 
obtained using the experimentally measured outer wall tem- 
peratures. These results are presented in tables III and VI. 
Table III contains the calculated nozzle inner wall tempera- 
tures, and table VI contains the calculated heat 11 ux to the nozzle 
wall. These temperatures and fluxes represent the experimen- 
tally determined values. The distribution of temperature along 
the length of the nozzle is shown in figure 28 for a typical firing 
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Propellant mixture ratio, O/F 


Figure 25.— Thrust chamber performance efficiency. Area ratio, e, 440:1 . 
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Vacuum thrust coefficient efficiency, ^ percent 



Propellant mixture ratio, 0/F 

Fiqure 26— Experimental and predicted nozzle vacuum thrust coefficient. Area ratio, 
e, 440:1. 
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Propellant mixture ratio, O/F 

Figure 27— Nozzle thrust coefficient efficiency. Area ratio, e, 440:1 . 
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(reading 577). Shown here are the measured outside wall 
temperatures and the calculated inner wall temperatures as a 
function of the nozzle length expressed as the expansion area 
ratio. Figure 29 shows the corresponding variation in calcu- 
lated heat flux. Although the heat flux decreases to compara- 
tively low values toward theexit of the nozzle, the corresponding 


nozzle surface areas involved increase substantially because of 
the nozzle contour. As a result, the heal rate to length ratio 
(product of heat flux times local circumference) of the trans- 
ferred heat nearer the exit becomes more significant than is 
apparent in figure 29, which considers only the heat flux vari- 
ation. This increased significance is apparent in figure 30, 
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Figure 28— Typical nozzle wall temperature distribution (reading 577). 
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Figure 29. — Typical calculated nozzle heat flux distribution (reading 577, uncorrected 
Q"). 
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Expansion area ratio, e 

Figure 30.— Typical nozzle heat rate/length distribution versus expansion area ratio 
(reading 577, uncorrected q"). 
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which is a plot of the heat rate to length ratio as a function of the 
axial length expressed in terms of the expansion area ratio for 
a typical firing (reading 577). In addition, the total heat rate 
from an area ratio of 140: 1 to an area ratio of 1025: 1 at the exit 
was calculated to show the total heat load for the nozzle skirt. 
This was done by integrating heat flux values over the wetted 
surface area of the nozzle (table VII). Figure 31 plots the heat 
rates as a function of the combustion chamber pressure. A 
faired line of the slope (P ,) 0 8 was drawn through the reading 
577 data point. For a given chamber configuration, heat- 
transfer rate can be considered proportional to the chamber 
pressure P raised to the 0.8 power (ref. 22). It is obvious that 
the data lie very parallel to this line, with a systematic scatter 
apparent because of the OIF variation of the firings. 

To reconcile the P variation of these data, the heat rate values 

< 

were adjusted to what they would have been if all the firings had 
been at the same P.. The P c selected was that of the typical 
firing, reading 577, which was 2063. 1 psia. The other heat rate 
values were corrected by multiplying them by (f\)° 8 . These 
values are tabulated in table VII and are also plotted in figure 32 
versus OIF . The data are well characterized by a straight line 
and show a minimal amount of scatter, which is caused pri- 
marily by experimental uncertainty. 

All these experimental results were then compared with the 
analytical predictions obtained from the TDK computer pro- 
gram as previously described. This computer code accounted 
for all the real effects expected with one exception. Because of 


the very specific nature of the various injectors and combustors 
used, the code was unable to account for combustion efficiency. 
Since we were able to measure this efficiency, an empirical 
correction was made to the heat flux calculations. This correc- 
tion took the experimental values of heat flux and increased 
them to what they would have been had we had 100-percent 
combustion efficiency. Characteristic exhaust velocity C* is 
proportional to the square root of the combustion gas tempera- 
ture; therefore, (r| c * )“ should vary directly with the combustion 
temperature and heal flux. Hence, 

<1\\m ~t ( 17 ) 

(%*)“ 

The empirical values of r^ * as a function of OIF were read 
from the best fit curve, equation (15) from figure 14. The 
adjustedexperimental heat flux values are tabulated in table XII. 
Table XII also includes heat flux values predicted by the TDK 
computer code with the turbulent BLM and M ABL modules. In 
addition, a laminar boundary layer module was run using the 
MABL module. For reference, results can be found in table 
XVII of appendix B. Because there was no significant differ- 
ence between the turbulent BLM and MABL modules, the 
remaining calculations and discussion are limited to results 
calculated with the MABL module. 



Figure 32. — Nozzle skirt heat rate (area ratio, e, 1 40:1 to 1025:1) as a function of 
mixture ratio (corrected to P Cj6 = 2063.1 psia). Curve represents heat rate, Q. 
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TABLE XII — EXPERIMENTAL <%-. = 100 PERCENT) AND TDK TURBULENT HEAT FLUX DISTRIBUTIONS 


Reading 

Effective 
combustion 
chamber total 
pressure at nozzle 
entrance. 

Propellant 

mixture 

ratio. 

O/F 

Expansion area ratio, £ 

50 

100 

200 | 

Heat flux to nozzle walls 

MPa 

psia 

kW/m 2 

2 

Btu/in."-s 

kW/tn 2 

2 

Blu/in."-s 

kW/m 2 

T 

Btu/in.‘-s 

Experimental adjusted to tv* = 100 percent 

569 

12.326 

1787.7 

3.89 

1428.00 

0.8738 

774.80 

0.4741 

413.46 

0.2530 

570 

12.645 

1834.0 

5.97 

1857.97 

1.1369 

1025.82 

.6277 

538.48 

.3295 

571 

12.488 

181 1.1 

4.70 

1643.72 

1 .0058 

902.76 

.5524 

470.83 

.2881 

575 

14.350 

2081.2 

4.65 

1815.81 

l.llll 

1009.31 

.6176 

547.80 

.3352 

576 

14.605 

2118.2 

5.68 

2013.55 

1.2321 

1117.66 

.6839 

598.62 

.3663 

577 

14.225 

2063. 1 

4.47 

1621.66 

.9923 

942.47 

.5767 

491.58 

.3008 

580 

16.364 

2373.3 

4.27 

1856.50 

1.1360 

1070.43 

.6550 

588.98 

.3604 

j TDK/MABL turbulent | 

569 

12.326 

1787.7 

3.89 

2556.94 

1 .5646 

1383.06 

0.8463 

714.66 

0.4373 

570 

12.645 

1834.0 

5.97 

3098.37 

1 .8959 

1670.53 

1 .0222 

853.89 

.5225 

571 

12.488 

1811.1 

4.70 

2754.85 

1 .6857 

1491.74 

.9128 

768.26 

.4701 

575 

14.350 

2081.2 

4.65 

3043.78 

1 .8625 

1673.79 

1.0242 

853.73 

.5224 

576 

14.605 

21 18.2 

5.68 

3372.76 

2.0638 

1833.13 

1.1217 

942.47 

.5767 

577 

14.225 

2063.1 

4.47 

2975.80 

1 .8209 

1618.23 

.9902 

829.05 

.5073 

580 

16.364 

2373.3 

4.27 

3242.83 

1 .9843 

1780.18 

1 .0893 

922.04 

.5642 


TDK/BLM turbulent j 

569 

12.326 

1787.7 

3.89 

2553.84 

1 .5627 

1361.98 

0.8334 

764.83 

0.4680 

570 

12.645 

1834.0 

5.97 ' 

3222.90 

1.9721 

1719.06 

1.0519 

956.69 

.5854 

571 

12.488 

1811. 1 

4.70 

2781.65 

1.7021 

1487.00 

.9099 

835.26 

.5111 

575 

14.350 

208 1 .2 

4.65 

3105.23 

1.9001 

1665.95 

1.0194 

926.45 

.5669 

576 

14.605 

21 18.2 

5.68 

3493.69 

2.1378 

1868.92 

1.1436 

1037.91 

.635 1 

577 

14.225 

2063. 1 

4.47 

301 1.75 

1.8429 

1603.85 

.9814 

892.30 

.5460 

580 

16.364 

2373.3 

4.27 

3245.94 

1 .9862 

1760.57 

1.0773 

984.80 

.6026 


Reading 

Expansion area ratio, £ 


300 


388 

500 

635 

800 

975 


Heat 11 ux to nozzle walls as measured 


kW/m 2 

2 

Btu/inT-s 

kW/m 2 

Btu/inT-s 

kW/m 2 

Btu/inT-s 

kW/m 2 

T 

Btu/inT-s 

kW/m 2 

Btu/in."-s 

kW/m 2 

Btu/in.^-s 




Experimental data adjusted to tv* = 100 percent 


569 

297.92 

0.1823 

239.42 

0.1465 

180.91 

0.1 107 

153.13 

0.0937 

143.49 

0.0878 

122.24 

0.0748 

570 

378.33 

.2315 

297.76 

.1822 

220.79 

.1351 

170.78 

.1045 

133.84 

.0819 

1 1 1 .62 

.0683 

571 

336.16 

.2057 

267.04 

.1634 

200.69 

1228 

155.74 

.0953 

126.33 

.0773 

103.77 

.0635 

575 

384.70 

.2354 

302.50 

.1851 

225.36 

.1379 

178.30 

.1091 

140.38 

.0859 

118.97 

.0728 

576 

414.12 

.2534 

327.67 

.2005 

243.50 

.1490 

189.90 

.1162 

145.77 

.0892 

124.20 

.0760 

577 

353.00 

.2160 

279.95 

.1713 

209.18 

.1280 

166.20 

.1017 

136.62 

.0836 

1 18.32 

.0724 

580 

411.67 

.2519 

327.34 

.2003 

245.79 

.1504 

189.25 

.1158 

148.55 

.0909 

127.47 

.0780 


TDK/MABL turbulent 1 

569 

469.03 

0.2870 

355.78 

0.2177 

263.77 

0.1614 

201.99 

0.1236 

150.84 

0.0923 

117.99 

0.0722 

570 

557.28 

.3410 

421.96 

.2582 

317.21 

.1941 

238.93 

.1462 

166.53 

.1019 

1 39.24 

.0852 

571 

505.47 

.3093 

381.92 

.2337 

288.28 

.1764 

216.54 

.1325 

162.61 

.0995 

126.98 

.0777 

575 

564.47 

.3454 

427.19 

.2614 

323.91 

.1982 

243.34 

.1489 

181.24 

1109 

141.85 

.0868 

576 

616.44 

.3772 

466.09 

.2852 

351.04 

.2148 

263.44 

.1612 

196.76 

.1204 

153.78 

.0941 

577 

548.94 

.3359 

414.61 

.2537 

313.28 

.1917 

236.64 

.1448 

176.50 

.1080 

138.42 

.0847 

580 

607.12 

.3715 

462.66 

.2831 

347.93 

.2129 

262.30 

.1605 

196.27 

.1201 

153.46 

.0939 



TDK/BLM turbulent 1 

569 

493.38 

0.3019 

360.51 

0.2206 

275.37 

0.1685 

206.90 

0.1266 

155.58 

0.0952 

122.08 

0.0747 

570 

615.13 

.3764 

449.58 

.2751 

341.23 

.2088 

257.39 

.1575 

193.00 

.1181 

151.82 

.0929 

571 

538.65 

.3296 

392.55 

.2402 

302.34 

.1850 

226.34 

.1385 

170.13 

.1041 

133.68 

.0818 

575 

601.08 

.3678 

437.00 

.2674 

334.04 

.2044 

251.51 

.1539 

188.59 

.1 154 

148.39 

.0908 

576 

667.43 

.4084 

488.64 

,2990 

372.28 

.2278 

278.80 

.1706 

209.5 1 

.1282 

164.24 

.1005 

577 

577.05 

.3531 

423.27 

.2590 

323.91 

.1982 

243.99 

.1493 

183.04 

.1120 

144.14 

.0882 

580 

642.91 

.3934 

470.66 

.2880 

358.72 

.2195 

269.16 

.1647 

202.48 

.1239 

158.69 

.0971 
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Figure 33 compares the TDK-predicted heat flux values with 
the commensurate experimentally measured heat flux values 
for reading 577 (adjusted forrj^*). This is a plot of the heat flux 
variation along the axial length of the nozzle expressed as the 
expansion area ratio. As seen in the figure, the experimental 
values fall below’ the heat flux values predicted by TDK for a 
turbulent boundary layer but are above the values predicted 
with a laminar boundary layer assumption. This was typical for 
all seven of the firings tabulated. This same variance between 
prediction and experiment is evident in figure 34, which is a plot 
of the heal rate to length ratio as a function of the location in the 
nozzle for reading 577. The heat rate to length ratio, which is the 
product of the heat flux times the local circumference, can take 
into consideration that the surface area of the nozzle is greater 
toward the exit because of the nozzle contour. Again, the 
experimental values fell below the heat rate to length ratios 
predicted by TDK for the turbulent boundary layer case but 
were above the values predicted with a laminar boundary layer 
assumption. The distributions of the heat flux and the heat rate 
to length ratio along the length of the nozzle were similar for all 
seven firings. Figure 35 illustrates the variation of the heat rate 
with respect to the OIF. This plot illustrates the total heat 
transferred from an area ratio of 140: 1 to 1025: 1 as a function 
of OIF. As with the previous two plots, the experimental results 
fall between the turbulent and laminar TDK predictions. Heat 


rate is the integral of the heat flux over the wetted surface area 
of the nozzle. All seven firings are represented by the data 
points obtained from table VII. Heat rates in this plot were 
generalized to a common combustion chamber pressure of 
2063. 1 psia, as they were in figure 32. The heat rates were also 
adjusted to the heat transfer values that would have occurred if 
the combustion were perfect (r| c * = I CX) percent) so that the 
TDK values could be compared to them. As in figures 33 
and 34, the experimentally obtained values fall below' the 
TDK-predicted values. This shortfall is quite consistent over 
the entire mixture ratio range, varying from 40 percent at 
OIF = 4 to 43 percent at OIF = 6. 

Boundary Layer 

Twelve firings were conducted at nominal combustion 
chamber pressures of 1800 psia, evenly divided amongst pro- 
pellant mixture ratios of 4, 5, and 6. Six of these firings were 
with the 1037: 1 -area-ratio nozzle, and six were with the nozzle 
truncated to the 440: 1 area ratio. In all these cases, boundary- 
layer pressure profile measurements were obtained with the 
previously described boundary-layer total pressure rakes. 

To verify that the thrust performance and combustion per- 
formance of the 12 firings were representative of the firings with- 
out the rakes, the C*, C FV , and / v values were calculated 



10 1 10 2 10 3 

Expansion area ratio, e 


Figure 33. — Typical calculated (experimental) and predicted nozzle heat flux distri- 
bution (reading 577). Experimental values are corrected for characteristic exhaust 
velocity efficiency. 
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Figure 34. — Typical calculated (experimental) and predicted nozzle heat rate/length 
distribution (reading 577). Values are corrected for characteristic exhaust velocity 
efficiency. 



Figure 35. — Calculated (experimental) and predicted nozzle skirt heat rate as a 
function of mixture ratio. Area ratio, e, 140:1 to 1025:1 . 
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and compared with the results of the firings without the rakes. 
Table XIII lists the thrust performance results of the 12 fir- 
ings, and figure 36 plots r| r * versus OIF. The faired curve 
represented by the equation shown is the curve fit for the 
data without the rakes previously described. The plotted data 
points all fall within the scatter from earlier firings. Figure 37 
shows the nozzle C f y data of table XIII plotted versus O/F. 
In spite of the drag produced by the boundary- layer rakes and 


the small difference in area ratio, the thrust performance was 
essentially the same as for the previous tests. This is obvious in 
comparing the thrust coefficient data points with the fitted 
curves. Agreement with the previous tests is also evident 
in figure 38, which is a plot of the T| r as a function of OIF. 
The faired lines represent the mean values of the previous 
firings. The data points, which are from table XIII, agree with 
the faired curves within the apparent scatter. 


TABLE XIII —RESULTS OF BOUNDARY LAYER TESTS 


Reading 

Nozzle throat area, 
A, 

Nozzle exit 
expansion 
area ratio. 

f 

Measured chamber pressure 

Effective combustion 
chamber 
total pressure at 

nozzle entrace,* 1 

Pc.c 

Propellant 
mixture ratio. 
O/F 

At injector end. 

Corrected tor momentum 
pressure loss. 

1 

cm" 

-> 

in." 

MPa 

psia 

MPa 

psia 

MPa 

psia 

589 

5.007 

0.776 

1037 

12.545 

1820.0 

12.507 

1814.4 

12.383 

1796.5 

3.97 

592 







13.008 

1887.1 

12.940 

1877.2 

12.787 

1855.1 

5.91 

593 







12.678 

1839.2 

12.621 

1831.0 

12.484 

1811.1 

4.98 

596 







12.731 

1847.0 

12.676 

1838.9 

12.539 

1819.1 

4.91 

597 







12.580 

1825! 

12.544 

1819.8 

12.421 

1802.0 

3.89 

598 







13.090 

1899.0 

13.022 

1889.1 

12.870 

1867.1 

5.80 

607 





440 

12.905 

1872.2 

12.836 

1862.2 

12.684 

1840.2 

5.99 

608 







12.651 

1835.3 

12.594 

1827.0 

12.458 

1807.3 

4.95 

609 







12.480 

1810.6 

12.443 

1805.1 

12.321 

1787.4 

3.91 

612 







12.978 

1882.8 

12.911 

1873.0 

12.760 

1851.1 

5.85 

613 







12.781 

1854.2 

12.727 

1846.4 

12.592 

1826.8 

4.76 

614 


r 

' 




12.576 

1824.5 

12.542 

1819.5 

12.420 

1801.8 

3.84 


Reading 

Vacuum thrust. 

E V 

Ambient pressure around 
nozzle. 

Pa 

Characteristic exhaust 
velocity, 

C* 

Characteristic 

exhaust 

velocity 

efficiency, 

lie- 

percent 

N 

ib, 

kPa 

psia 

m/s 

ft/s 

589 

1 1 885 

2671.9 

1.313 

0.1905 

2497 

8193 

99.9 

592 

12 920 

2904.6 

1.502 

.218 

2341 

7680 

98.9 

593 

12 318 

2769.4 

1 .378 

.200 

2418 

7934 

99.2 

596 

12 369 

2780.9 

1 .450 

.210 

2420 

7938 

99.1 

597 

1 l 860 

2666.4 

1 .4(H) 

.203 

2492 

8175 

99.5 

598 

13 015 

2926.1 

1 .540 

.223 

2349 

7707 

98.9 

607 

12 278 

2760.4 

.892 

.129 

2350 

7710 

99.6 

608 

1 1 777 

2647.6 

.807 

.117 

2448 

8031 

1 00.4 

609 

1 1 298 

2540.0 

1.668 

.242 

2514 

8248 

100.5 

612 

12 427 

2793.8 

.928 

.135 

2372 

7782 

100.0 

613 

11 810 

2655.1 

.885 

.128 

2458 

8065 

100.2 

614 

1 1 322 

2545.5 

.849 

.123 

2509 

8231 

100.1 


Calculated with low nozzle exit expansion area ratio correlation. 
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TABLH XIII. — Concluded. 


Reading 

Fuel injection 

Oxidizer injection 

Propellant flow rate, 
rii 

Pressure, 

Pr, 

Temperature, 

Tfi 

Pressure. 

Pr, 

Temperature, 

Tf, 


MPa 

psia 

K 

°R 

MPa 

psia 

K 

°R 

kg/s 

lb„/s 

589 

16.413 

2381.1 

299.4 

539.0 

13.670 

1983.2 

111.8 

201.3 

2.48 

5.48 

592 

15.479 

2245.6 

300.7 

541.3 

14.490 

2102.2 

109.3 

196.7 

2.74 

6.03 

593 

15.619 

2265.9 

300.1 

540.1 

13.962 

2025.6 

109.1 

196.3 

2.59 

5.70 

596 

15.738 

2283.2 

299.3 

538.8 

14.044 

2037.4 

118.3 

212.9 

2.60 

5.72 

597 

16.569 

2403.8 

298.2 

536.8 

13.678 

1984.3 

1 13.7 

204.6 

2.50 

5.50 

598 

15.637 

2268.5 

298.8 

537.9 

14.607 

2119.1 

107.1 

192.7 

2.74 

6.05 

607 

15.245 

221 1.6 

298.3 

536.9 

14.367 

2084.3 

112.4 

202.4 

2.70 

5.96 

608 

15.510 

2250. 1 

298.2 

536.7 

13.902 

2016.9 

1 10.9 

199.7 

2.55 

5.62 

609 

16.303 

2365.1 

298.0 

536.4 

13.549 

1965.6 

1 1 1 .9 

201.5 

ir, 

-t 

ri 

5.41 

612 

15.399 

2234.0 

301.2 

542.1 

14.449 

2096.2 

lll.l 

199.9 

2.69 

5.94 

613 

15.847 

2299.0 

300.3 

540.6 

14.042 

2037.2 

1 10.9 

199.6 

2.57 

5.66 

614 

16.582 

2405.6 

299.9 

539.8 

1 3.648 

1980.0 

1 14.2 

205.6 

2.48 

5.47 


Reading 

Measured 
vacuum thrust 
coefficient. 

Q.v 

Vacuum 

thrust 

coefficient 

efficiency, 

V • 
vp.v 

percent 

Vacuum 

specific 

impulse. 

!sp,V» 

s 

Vacuum 

specific 

impulse 

efficiency, 

n, 

Up.v 

percent 

589 

1.917 

97.4 

488.0 

97.4 

592 

2.018 

96.3 

481.6 

95.3 

593 

1.971 

97.0 

485.9 

96.2 

596 

1.970 

97.2 

486.0 

96.2 

597 

1 .907 

97.2 

484.5 

96.8 

598 

2.020 

96.8 

483.8 

95.7 

607 

1 .933 

93.4 

463.2 

93.0 

608 

1.888 

94.1 

471.2 

94.5 

609 

1.831 

94.2 

469.5 

94.6 

612 

1 .945 

94.3 

470.5 

94.4 

613 

1 .873 

93.9 

469.5 

94.1 

614 

1.821 

93.8 

465.7 

93.9 



Figure 36. — Characteristic exhaust velocity efficiency as a function of mixture ratio for 
tests with boundary layer rakes installed. 
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Vacuum thrust coefficient efficiency, , percent 



3 4 5 6 7 

Propellant mixture ratio, 0/F 

Figure 37. — Nozzle thrust performance as a function of mixture ratio for tests with 
the boundary layer rakes installed. 



3 4 5 6 7 

Propellant mixture ratio, O/F 

Figure 38. — Nozzle thrust coefficient efficiency as a function of mixture ratio for tests 
with the boundary layer rakes installed. 


38 


NASA/TP — 1999-208522 




Figure 39.— Overall thruster efficiency as a function of mixture ratio for tests with the 
boundary layer rakes installed. 



Figure 40. — Overall thruster performance as a function of mixture ratio for tests with 
the boundary layer rakes installed. 
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Overall performance of the thruster is illustrated by figure 39. 
which is a plot of r} 7 ^ v ,as a function of OIF. Faired curves of 
the mean values of die previous firings are shown with the data 
points of table XIII superimposed; again, there is very good 
agreement within the scatter band. Figure 40 shows the actual 
specific impulse attained versus OIF. Here, the faired curve 
shows the mean values for the previous nonrake firings, and the 
individual data points are the values for firings with the rakes. 

The total pressure measurements, along with pertinent test 
conditions, are listed in table XIV, which shows test results for 
both the 1037: 1- and 440:1 -area-ratio nozzles. Figure 41 shows 
a graphic illustration of the pressure profile measured for the 


1037:1 configuration, along with bands of boundary layer 
thickness 5 and displacement thickness 8* from the TDK 
predictions. To avoid the scatter caused by slight variances in 
combustion chamber pressure among the six firings plotted, the 
measured total pressures were normalized by the combustion 
chamber pressure P ce and the resultant pressure ratio was 
plotted versus the distance from the wall. The symbols for the 
data points are coded to indicate the nominal OIF l or each point. 
A slight dependence on mixture ratio is evident in the distribu- 
tion on the plot. Also evident is the consistent agreement of the 
data: all the data obtained in six separate firings with three 
different total pressure rakes defined the same total pressure 


TABLF XIV.— EXPERIMENTALLY MEASURED BOUNDARY LAYER TOTAL PRESSURE 


Reading 

Nozzle exit 
expansion area 
ratio. 

F 

Effective combustion 
chamber total pressure at nozzle 

a 

entrance. 

Pee 

Propellant mixture 
ratio, 

O/E 

Rake total pressure 

Tube A ( 1 in.) 

Tube A (4 in.) 

Exact area ratio of tube. 1019.6 
(0.0440 in. from nozzle wall) 

Exact area ratio of tube. 1 0 1 9.7 
(0.0655 in. from nozzle wall) 

MPa 

psia 

kPa 

psia 

kPa 

psia 

589 

1037 

12.383 

1796.5 

3.97 



9.263 

1.344 

592 



12.787 

1855.1 

5,91 



8.546 

1.240 

593 



12.484 

1811.1 

4.98 



8.403 

1.219 

596 



12.539 

1819.1 

4.91 

5.714 

0.829 



597 



12.421 

1802.0 

3.89 

5.985 

.868 



598 


f 

12.870 

1867.1 

5.80 

5.692 

.826 









Exact area ratio of tube. 418.7 

Exact area ratio of tube. 41 8.9 







(0.0680 in. from nozzle wall ) 

(0,0900 in from nozzle wall) 







kPa 

psia 

kPa 

psia 

607 

440 

12.684 

1840.2 

5.99 



48.499 

7.036 

608 



12.458 

1807.3 

4.95 



43.385 

6.294 

609 



12.321 

1787.4 

3.91 



38.842 

5.635 

612 



12.760 

185 LI 

5.85 

44.625 

6.474 



613 



12.592 

1826.8 

4.76 

43.564 

6.320 



614 


f 

12.420 

1801.8 

3.84 

43.154 

6.261 




Reading 




Rake total 

pressure 



1 


Tube B ( 1 in.) 

Tube C ( 1 in.) 

Tube D ( 1 in.) 

Tube B (4 in.) 


Exact area ratio of tube. 1019.6 
(0.0895 in. from nozzle wall) 

Exact area ratio of tube. 1020.5 
(0.3845 in. from nozzle wall) 

Exact area ratio of tube, 1020.7 
(0.4380 in. from nozzle wall) 

Exact area ratio of tube, 1021 .3 
(0.6796 in. from nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

589 







22.177 

3.217 

592 







20.969 

3.042 

593 







20.758 

3.011 

596 

11.084 

1.608 

17.653 

2.561 

19.575 

2.840 



597 

1 1.429 

1.658 

18.073 

2.622 

20.221 

2.934 



598 

11.091 

1.609 

18.122 

2.629 

19.679 

2.855 




Exact area ratio of tube. 420. 1 
(0.2845 in. from nozzle wall) 

Exact area ratio of tube, 42 1.0 
(0.4195 in. from nozzle wall) 

Exact area ratio of tube, 422.5 
(0.6440 in. from nozzle wall) 

Exact area ratio of tube, 422.9 
(0.7045 in from nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

607 







49.457 

7.175 

608 







48.765 

7.075 

609 







47.098 

6.833 

612 

66.550 

9.655 

78.547 

1 1 .395 

88.351 

12.818 



613 

65.898 

9.560 

81.269 

1 1 .790 

88.403 

12.825 



614 

64.621 

9.375 

79.177 

11.487 

88.342 

12.816 




'Damaged tube. 
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TABLED XIV — Concluded. 


Reading 

Rake total pressure 


Tube E ( 1 in.) 

Tube E ( 1 in.) 

Tube G (1 in. ) 

lube C (4 in.) 


Exact area ratio of tube. 1021.5 

Exact area ratio of tube. 1 02 1 .7 

Exact area ratio of tube. 1022.4 

Exact area ratio of tube. 1022.9 


(0.7535 in. from nozzle wall ) 

(0.8240 in. from nozzle wall ) 

(1.1 280 in. from nozzle wall) 

( 1.3135 in fri 

>m nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 



S89 




'a ifl-fifi 11 '- if 





Sep 

















4.579 

596 

24.148 

3.503 

25.719 


30.490 

4.423 



597 

25.062 

3.636 

26.863 


31.898 

4.628 



598 

24.565 

3.564 

22.626 


30.687 

4.452 




Exact area rat it 

> of tube. 423. 1 

Exact area ratio of tube. 424.9 

Exact area ratio of tube. 425.3 

Exact area ratio ot tube. 426.8 


(0.7455 in. from nozzle wall) 

( 1 .0125 in. from nozzle wall ) 

( 1 .0805 in. from nozzle wall) 

( 1 .2970 in. fr< 

>m nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

607 







61.990 

8.993 

608 







59.774 

8.672 

609 







58.040 

8.420 

612 

88.322 

12.813 

52.079 

;, 7.555 

70.708 

a 10.258 



6 1 

55 1 

7 997 

25 7 1 9 

A ^ 7^| 

30.488 

a 4 4~>3 



614 

56.512 

8.199 

51.215 

‘*7.430 

39. 1 38 

‘'5.678 




Reading 

Rake total pressure 


Tube L) (4 in. ) 

lube E (4 in.) 

l ube E (4 in. ) 

l ube G (4 in.) 


Exact area ratio of tube. 1024.5 

Exact area ratio of tube, 1026. 1 

Exact area ratio of tube, 1027.7 

Exact area ratii 

» of tube. 1029 3 


( 1.9345 in. from nozzle wall ) 

(2.5645 in. from nozzle wall) 

(3.1615 in. from nozzle wall) 

(3.7865 in. fn 

>ni nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

589 

38. 169 

5.537 

32.91 1 

4.775 

31.239 

4.532 

28.477 

4 131 

592 

36.235 

5.257 

30.38 1 

4.408 

31 .069 

4.507 

28.727 

4.168 

593 

36.281 

5.263 

31.744 

4.605 

30.946 

4.490 

28.290 

4.104 

596 









597 

598 










Elxact area ratic 

of tube. 430.8 

Exact area rat 

io of tube. 435.0 

Exact area ratio of tube. 439. 1 

Exact area ratio of tube. 443.2 1 


( 1 .9155 in. from nozzle wall) 

(2.5375 in. from nozzle wall) 

(3. 1 5 1 5 in from nozzle wall) 

(3.7625 in trc 

>m nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

607 

74.354 

10.787 

82.020 

1 1 .899 

88.166 

12.791 

45.489 

*6.599 

608 

77.799 

1 1.287 

83.583 

12.126 

88.192 

12.795 

44.363 

'6.435 

609 

75.397 

10.938 

85.257 

12.369 

88.384 

12.822 

43.1 12 

*6.254 

6 1 2 









6 1 3 
614 










* Damaged lube. 


profile. Figure 42 illustrates the pressure profile obtained from 
the 440: 1 configuration tests and shows hands of 5 and 5* from 
the TDK predictions. A slight dependence on mixture ratio is 
alsoevident here as is the consistent agreement of the six firings 
and three rakes, defining one total pressure profile at the 440: 1 
location. 

For comparison, theoretical predictions were obtained from 
the TDK computer axle for the 12 firings. A program option 
was used that provides the results that would be created down- 
stream of a normal shock if a pitot tube was placed into the How 
field. Table XV tabulates these results, and they are graphically 
displayed as solid lines in figures 41 and 42 for the 1037:1 -and 
440: 1 -expansion-area-ratio nozzles, respectively. As figures 4 1 


and 42 illustrate, there is good agreement between the experi- 
mental results and the analytical predictions for the first 
2.54 cm (1.0 in.) and 1.27 cm (0.5 in.) in the 1037:1 and 440:1 
nozzles, respectively. However, the TDK code overpredicted 
the measured free-stream pressure at both axial locations, 
perhaps because of combustion or shock losses in the nozzle. 
The free-stream pressure in these exit planes is a function of the 
radius, making it difficult to sort out experimentally the free- 
stream viscous pressure loss from the shock and combustion 
losses and making it difficult to measure the edge of the 
boundary layer. It is impossible to translate the measured 
global combustion losses into the local-exit-plane pressure 
losses. The experimental results indicate very little variation 
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TABLE XV — TDK/MABL PREDICTIONS OF BOUNDARY LAYER TOTAL PRESSURE 


Reading 

Nozzle exit 
expansion 
area ratio, 
£ 

Effective combustion 
chamber total pressure 

at nozzle entrance, 11 

Pc.e 

Propel la 
mixtur 
ratio, 
O/F 

MPa 

psia 

589 

1037 

12.383 

17%.5 

3.97 

592 



12.787 

1855.1 

5.91 

593 



12.484 

1811.1 

4.98 

596 



12.539 

1819.1 

4.91 

597 



12.421 

1802.0 

3.89 

598 


F 

12.870 

1867.1 

5.80 


607 

440 

12.684 

1 840.2 

5.99 

608 



12.458 

1807.3 

4.95 

609 



12.321 

1787.4 

3.91 

612 



12.760 

1851.1 

5.85 

613 



12.592 

1826.8 

4.76 

614 



12.420 

1801.8 

3.84 


TDK/M ABL rake total pressure predictions 


Tube A ( 1 in.) 


Exact area ratio of tube, 1019.6 
(0.0440 in. from nozzle wall) 


kPa 


5.185 

4.853 

5.347 


psia 


0.752 

.704 

.776 


Exact area ratio of lube, 4 1 8.7 
(0.0680 in. from nozzle wall) 


kPa 


psia 


Tube A (4 in. ) 


Exact area ratio of tube, 1019.7 
(0.0655 in. from nozzle wall) 


kPa 


6.757 
7.631 
7 118 


psia 


0.980 
1.107 
1 .033 


Exact area ratio of tube. 418.9 
(0.0900 in from nozzle wall) 


kPa 


psia 


43.054 

40.689 

38.835 


6.246 

5.903 

5.634 


48.051 

45.514 

44.391 


6.971 

6.603 

6.440 


Reading 


589 

592 

593 

596 

597 

598 


TDK/MABL rake total pressure predictions 


Tube B ( 1 in. ) 


Exact area ratio of tube. 1019.6 
(0.0895 in. from nozzle wall) 


kPa 


9.002 

8.561 

9.308 


psia 


1 306 
1 .242 
1 .350 


Exact area ratio of tube. 420. 1 
(0.2845 in. from nozzle wall) 


kPa 


psia 


Tube C ( 1 in.) 


Exact area ratio of tube, 1()20.5| 
(0.3845 in. from nozzle wall 


kPa 


18.935 

18.500 

19.414 


psia 


2.747 

2.684 

2.817 


Exact area ratio of tube, 42 1 .0 
(0.4195 in. from nozzle wall ) 


kPa 


psia 


Tube D ( I in.) 


Exact area ratio of tube. 1020.7 
(0.4380 in. from nozzle wall) 


kPa 


19.933 

19.583 

20.371 


psia 


2.892 

2.841 

2.955 


Exact area ratio of tube, 422.5 
(0.6440 in. from nozzle wall) 


kPa 


psia 


Tube B (4 in.) 


Exact area ratio of tube, 1 02 1 .3 
(0.6796 in from nozzle wall) 


kPa 


24.346 

24.525 

24.098 


psia 


3.532 

3.558 

3.496 


Exact area ratio of tube. 422.9 
(0.7045 in. from nozzle wall) 


kPa 


psia 


607 

608 
609 
612 

613 

614 


70.922 

70.129 

69.895 


10.289 

10.174 

10.140 


82.826 

83.536 

82.578 


12.016 
12.119 
I 1 .980 


98.742 

102.313 

106.842 


14.325 

14.843 

15.500 


101 .437 
104.643 
111.232 


14.716 

15.181 

16.137 


1 Outside of theoretical boundary layer. 
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TABLE XV — Concluded. 


Reading 

TDK/M ABL rake total pressure predictions 


Tube E ( 1 in.) 

Tube F ( 1 in.) 

Tube G ( i in.) 

Tube C (4 in.) 


Exact area ratio of tube, 102 1 .5 
(0.7535 in. from nozzle wall) 

Exact area ratio of tube. 1 02 1 .7 
(0.8240 in. from nozzle wall) 

Exact area ratio of tube, 1022.4 
( 1.1280 in. from nozzle wall) 

Exact area ratio of tube. 1022.9 
(1.31 35 in. from nozzle wall ) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

589 







36.650 

5.317 

592 







34.1 13 

4.949 

593 







34.688 

5.032 

596 

25.621 

3.717 

26.862 

3.897 

32.025 

4.646 



597 

25.980 

3.769 

27.413 

3.977 

33.403 

4.846 



598 

25.704 

3.729 

26.848 

3.895 

31.570 

4.580 




Exact area ratio of tube, 423. 1 
(0.7455 in. from nozzle wall) 

Exact area ratio of tube, 424.9 
( 1 .0125 in. from nozzle wall) 

Exact area ratio of tube, 425.3 
( 1.0805 in. from nozzle wall) 

Exact area ratio of tube, 426.8 
( 1.2970 in. from nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

607 







61.990 

8.993 

608 







59.774 

8.672 

609 







58.040 

8.420 

612 

HU. 774 

15.200 

(a) 

(a) 

(a) 

(a) 



613 

109.082 

15.825 

(a) 

(a) 

(a) 

(a) 



614 

1 15.044 

16.690 

(a) 

(a) 

(a) 

(a) 




Reading 

TDK/M ABL rake total pressure predictions 


Tube D (4 in.) 

Tube E (4 in.) 

Tube F (4 in.) 

Tube G (4 in.) 


Exact area ratio of tube, 1024.5 
{ 1 .9345 in. from nozzle wall) 

Exact area ratio of tube, 1 026. 1 
(2.5645 in. from nozzle wall) 

Exact area ratio of tube, 1027.7 
(3.1615 in. from nozzle wall) 

Exact area ratio of tube. 1029.3 
(3.7865 in. from nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

589 

47.624 

6.909 

49.457 

7.175 

49.361 

7.161 

(a) 

(a) 

592 

42.027 

6.097 

43.771 

6.350 

43.660 

6.334 

(a) 

(a) 

593 

43.908 

6.370 

45.742 

6.636 

45.701 

6.630 

(a) 

(a) 

596 









597 









598 










Exact area ratio of tube, 430.8 
( 1.9155 in. from nozzle wall) 

Exact area ratio of tube, 435.0 
(2.5375 in. from nozzle wall) 

Exact area ratio of tube. 439. 1 
(3.1515 in. from nozzle wall) 

Exact area ratio of tube, 443.2 
(3.7625 in. from nozzle wall) 


kPa 

psia 

kPa 

psia 

kPa 

psia 

kPa 

psia 

607 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

608 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

609 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

612 









613 









614 










J Outside of theoretical boundary layer. 


due lo mixture ratio; they show a sudden drop in pressure and 
then a gradual asymptote to some core How value. The analyti- 
cal predictions, however, indicate a larger variation in the total 
pressure profile due to mixture ratio and also indicate a seem- 
ingly sharp knee where they asymptote to some significantly 
higher core flow value. The code seems to overpredict the 
boundary layer thickness. 

Within measurement accuracy, these data suggest a turbu- 
lent boundary layer profile. However, a subtle profile change 
developed between the 440:1 position and 1037:1, indicating 
that the near- wall profile was becoming more laminarlike, 


possibly in a relaminarization process caused by the highly 
favorable pressure gradient or, perhaps, as a laminar sublayer. 
This process, which the code did not predict, would account for 
the experimental performance and heat-transfer data falling 
between laminar and turbulent predictions. Predicting the meas- 
ure of laminarization then becomes important for predict- 
ing performance where laminar and turbulent predictions vary 
3 percent, and it becomes critical to predicting heat transfer 
where experimental measurements and laminar and turbulent 
predictions vary approximately 1 00 percent. 
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Concluding Remarks 

The results show that, based on the assumption oflurbulent 
How due to the throat Reynolds number, the TDK code predicts 
performance better for the 440: 1 -nozzle configuration. Refer- 
ence 6 results are considered to have val idaled the TDK/M ABL 
code up to an area ratio of 300: 1 . Therefore, on the basis of the 
current results and previous work, it would be reasonable to 
conclude that the TDK/MABL code is validated up to an area 
ratio of 440:1. The results from the TDK turbulent model 
provide a nearly 2- percent overprediction of experimental 
results at an area ratio of 1025:1 and only an approximately 
0.5-percent underprediction from the TDK laminar model. 
However, the TDK turbulent model underpredicted perfor- 
mance by only 0.25 percent for an area ratio o! 440: 1 . whereas 
the TDK laminar prediction was nearly 2.0-percent higher than 
experimental results. Although a 2-percent difference between 
the turbulent prediction and experiment is not desirable at the 
1025: 1 area ratio, it is better to have a code that underpredicts 
performance rather than one that is overly optimistic. Because 
the code seems to underpredicl performance at high area ratios, 
it is possible that the boundary layer growth is overstated. As 
the boundary layer grows, it reduces the in viscid core size, 
which might result in lower-than-expected predicted perfor- 
mance. The exact nature of the boundary layer flow phenom- 
ena is still one of the least understood and most difficult 
portions to model of any nozzle flow. Variations in the boundary 
model have generally been thought of as small for performance 
calculations; but for high-area-ratio nozzles, that is not the case. 

Summary of Results 

Experimentally attained rocket performance was compared 
with the current JANNAF methodology of performance pre- 
diction. A gaseous hydrogen and liquid oxygen rocket thruster 
with a 1 -in. -diameter throat was test fired at altitude, and the 
thrust performance, heat transfer rate, and total pressure pro- 
files were measured. Firings with both a 1025:1 and a truncated 
440: 1 -expansion-area-ratio configuration were compared 
with the predictions. The test firings were at combustion 
chamber pressures of 1800 to 2400 psia. and at propellant 
mixture ratios of 3.9 to 6.0. The Reynolds number, based on 
throat diameter, of the How through the throat was 1 .43xl0 6 
to 2.05x1 0 6 , depending on the mixture ratio and combustion 
chamber pressure. 

Two performance predictions were made with the TDK com- 
puter code for each of the firings, one each with the BLM 
boundary layer module and the M ABL boundary layer module. 
Differences between the turbulent BLM and MABL results 
proved insignificant, and further discussion was limited to the 
MABL results. Results were also compared with predictions 
using the laminar boundary layer model in the MABL module. 
Comparisons of predicted performance to experimentally 
attained thrust performance indicated that the experimentally 


attained performance was approximately 2.0-percent higher 
than the turbulent prediction and approximately 0.5-percent 
lower than the laminar prediction for the 1025:1 configura- 
tions. However, for the 440:1 configuration, the experimen- 
tally attained performance was approximately 0.25-percent 
higher than the turbulent prediction and approximately 
2.0-percent lower than the laminar prediction. 

Nozzle wall temperatures were measured on the outside of 
a thin-walled heat sink nozzle during the test firings. Nozzle 
heat fluxes were calculated from the time histories of these 
temperatures. Values of the integral of heat flux as a function 
of nozzle surface area were also calculated. So that they could 
be directly compared with the analytical predictions, the 
experimental values were adjusted in a number of ways. The 
heat flux, heat rate per unit length, and heat rale values were 
adjusted to what they would have been with complete combus- 
tion by the square of the characteristic exhaust velocity effi- 
ciency (n c D 2 - The heat rate values were adjusted to a uniform 
combustion chamber pressure P ( of 2063. 1 psia (reading 577) 
by the factor (P.) a8 . As a result, two systematic causes of data 
scatter were reconciled, and true variations in heat transfer as 
a result of other functions became obvious. A comparison of 
experimental heat rate to the analytical predicted values shows 
a very similar OIF dependence, although the experimental 
values are lower than the predicted values. This shortfall of 
experimentally measured heat transfer is also evident in the 
comparison of the experiment to the analysis of heal flux and 
heat rate per axial length. 

A separate series of high-pressure rocket nozzle firings were 
conducted to document the boundary layer profile of a high- 
area-ratio nozzle. The nozzles had expansion area-ratios of 
1037: 1 and 440: 1 with a nominal throat diameter of 2.54 cm 
( 1 .0 in.). Characteristic exhaust velocity, nozzle thrust coeffi- 
cient, and thruster specific impulse were determined and com- 
pared with nearly identical firings without boundary layer 
rakes to ensure applicability. As indicated by a comparison ol 
boundary layer total pressure profiles and analytical predic- 
tions. there was good agreement for 0.5 in. from the nozzle 
wall; but the further into the core that flow measurements were 
taken, the more TDK overpredicted the boundary layer thick- 
ness. Several possible explanations exist, such as the possibil- 
ity of relaminarization due to the highly favorable pressure 
gradient, the size of the laminar sublayer, or the type of 
turbulence waves present; however, more investigation is 
required. The difference between measured and predicted 
freestream pressures also indicates that local flow properties 
are significantly affected by combustion efficiency and shock 
losses. The current methodology does not account for this 
process, which may be a key to improving high-area-ratio 
performance predictions. 

Glenn Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio. March 31, 1999 
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Appendix A 
Symbols 


A , 

\ 


nozzle exit area, m 2 (in. 2 ) 
nozzle surface area, m 2 (in. 2 ) 
nozzle throat area, m 2 (in. 2 ) 
venturi throat area, m 2 (in. 2 ) 


k conductivity, W/m-K (Btu/sec-in.-°F) 

Lj starting axial position, m (in.) 

ending axial position, m (in.) 
m propellant mass flow 


C # 


C 


IV 


venturi discharge coefficient, dimensionless 
vacuum thrust coefficient, dimensionless 


O/F 


C i v 77 }{ODE) theoretical, one-dimensional-equilibrium 
(ODE) vacuum thrust coefficient (obtained 
from the Chemical Equilibrium Composition 
(CEC) program), dimensionless 

Cf v Th(TDK) theoretical, two-dimensional-kinetics (TDK) 
vacuum thrust coefficient, dimensionless 


C* 


C*monE) 


F 

F. 




characteristic exhaust velocity, m/s (fl/s) 

theoretical, one-dimensional -equilibrium 
characteristic exhaust velocity (obtained from 
the CEC program), m/s (ft/s) 

diameter, m (in.) 

thrust (corrected for aneroid effect), N (lb t ) 

vacuum thrust (experimentally measured thrust 
corrected to vacuum conditions), N (lb f ) 

acceleration due to gravity, 9.807 m/s 2 
(32.174 ft/s 2 ) 

proportionality constant, 1 kg-m/N-s 2 
(32.2 Ib m -ft/lb r s 2 ) 

theoretical subsonic specific impulse inside 
combustion chamber (obtained from the CEC 
program), N-s/kg (lb r s/lb m ) 


P,e 


P fi 


P/Pt 

Q 

<//' 

Re a 


propellant mixture ratio (oxidizer flow divided 
by fuel flow), dimensionless 

ambient pressure in test capsule, kPa (psia) 

chamber pressure 

static pressure at injector end of combustion 
chamber, MPa (psia) 

effective combustion chamber total pressure at 
nozzle entrance, MPa (psia) 

combustion chamber total pressure after com- 
bustion (P corrected for momentum pressure 
loss), MPa (psia) 

fuel injection pressure, MPa (psia) 

oxidizer injection pressure. MPa (psia) 

static pressure in nozzle, kPa (psia) 

static-to-total pressure ratio in combustion 
chamber (obtained from the CEC program), 
dimensionless 

heat rate, W (Btu/sec) 

heat flux to inner wall of nozzle, W/m 2 
(Btu/in. 2 -sec) 


Reynolds number based on momentum 
thickness 


I 


sp. V 


vacuum specific impulse, N-s/kg (lb^s/lb m ) R. 


local nozzle inner wall radius, m (in.) 


I S j } v jfiioDi) theoretical, one-dimensional equilibrium vac- 
uum specific impulse (obtained from the 
CEC program), N-s/kg (lb r s/lb m ) T 


local nozzle outer wall radius, m (in.) 
temperature, K (°R) 
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fuel injection temperature. K (°R) 


8 * 


boundary layer displacement thickness, cm (in.) 



a 

AP 

8 



velocity through venturi throat, m/s (in./s) 

mass-averaged injection velocity of propellants, 
m/s (ft/s) 

diffusivity, rrr/s 

nominal pressure drop, kPa (psid) 
boundary layer thickness, cm (in.) 


’iQ-.v 

%,,.v 

P 

c 

0 


nozzle exit expansion area ratio, A e JA r 
dimensionless 

thruster contraction area ratio, dimensionless 

characteristic exhaust velocity efficiency, 
percent 

vacuum thrust coefficient efficiency, percent 
vacuum specific impulse efficiency, percent 
fiuid density, kg/nr 1 (lb m /in.^) 
standard deviation, dimensionless 
angle between nozzle wall and axis, deg 
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Appendix B 

Laminar Boundary Layer Results 


Table XVI presents the TDK/MABL predictions with the laminar boundary layer condition, and table XVII compares 
the laminar boundary layer heat flux results with the TDK/MABL predictions. 


TABLE XVI.— TDK/MABL LAMINAR PREDICTIONS 


Reading 

Nozzle exit 
expansion 
area ratio, 
£ 

Effective combustion 
chamber total pressure at 
nozzle entrance. 

Propellant 
mixture ratio, 
O/F 

Predicted propellant flow rate 

MPa 

psia 

kg/s 

lb„/s 

569 

1025 

12.326 

1787.7 

3.89 

2.4926 

5.4961 

570 



12.645 

1834.0 

5.97 

2.7093 

5.9740 

571 



12.488 

181 LI 

4.70 

2.5718 

5.6709 

575 



14.350 

2081.2 

4.65 

2.95 1 1 

6.507 1 

576 



14.605 

21 18.2 

5.68 

3.0960 

6.8266 

577 



14.225 

2063.1 

4.47 

2.91 14 

6.4196 

580 



16.364 

2373.3 

4.27 

3.3321 

7.3472 

601 

440 

12.768 

1851.8 

6.15 

2.7196 

5.9968 

602 



12.542 

1819.0 

5.1 1 

2.5830 

5.6955 

603 



12.457 

1806.7 

4.01 

2.4944 

5.5001 


Reading 

Nozzle exit 

Computer code 


expansion 
area ratio. 

TDK/MABL Laminar 


£ 

Predicted characteristic 

Predicted vacuum thrust. 

Predicted 

Predicted 




exhaust velocity. 

Fv 

vacuum thrust 

vacuum thrust 




c* 

coefficient. 

coefficient 








Ci > 

efficiency. 




m/s 

ft/s 

N 

lb t 


V 

kj .v 
percent 

569 

1025 

2505.14 

8219.37 

11975 

2692.29 

1.92 

97.73 

570 



2364.44 

7757.74 

13088 

2942.47 

2.04 

97.32 

571 



2459.70 

8070.27 

12432 

2795.05 

1.96 

97.61 

575 



2463.32 

8082.16 

14271 

3208.30 

1.96 

97.75 

576 



2389.69 

7840.57 

14995 

3371.15 

2.03 

97.56 

577 



2475.18 

8121.07 

14070 

3163.33 

1.95 

97.76 

580 



2487.78 

8162.41 

16094 

3618.23 

1.94 

97.84 

601 

440 

2349.84 

7709.83 

12733 

2862.63 

1.99 

95.77 

602 



2430.36 

7974.01 

12137 

2728.56 

1 .93 

95.93 

603 



2499.57 

8201.08 

1 1670 

2623.55 

1.87 

95.96 
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TABLE XVI — Concluded 


Reading 

Nozzle exit 
expanison 
area ratio, 
£ 

Computer code 

ODL 

ODK 

MOC 

TDK/MABL Laminar 

Predicted 
vacuum specific 
impulse, 

Isp.V' 

sec 

Predicted 

vacuum 

specific 

impulse 

(adjusted). 

^sp.V’ 

sec 

Predicted 
vacuum specific 
impulse 
efficiency 
(adjusted). 

Lp.v- 

percent 

569 

1025 

500.63 

499.61 

494.93 

489.86 

488.97 

97.67 

570 



505.53 

502.12 

498.62 

492.55 

487.97 

96.53 

571 



504.43 

502.71 

498.44 

492.88 

491.16 

97.37 

575 



503.90 

502.51 

498.18 

493.05 

491.41 

97.52 

576 



505.55 

503.09 

499.37 

493.83 

490.05 

96.94 

577 



503.48 

502.22 

497.79 

492.75 

491.35 

97.59 

580 


r 

502.54 

501.62 

497.07 

492.46 

491.29 

97.76 

601 

440 

497.81 

494.36 

482.12 

477.36 

472.40 

94.90 

602 



498.99 

496.98 

483.61 

479.07 

476.72 

95.54 

603 



496.58 

495.52 

481.11 

477.00 

476.06 

95.87 
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TABLH XVII. — EXPERIMENTAL (iy. = 100 PERCENT) AND TDK/MABL LAMINAR HEAT FLUX DISTRIBUTIONS 


Reading 

Effective combustion 
chamber total pressure 
at nozzle entrance, 

Pc.e 

Propellant 
mixture ratio. 
O/F 

Expansion area ratio, e 

S» 

100 

200 

Heat flux to nozzle wall 

MPa 

psia 

kW/m 

Btu/in. '-s 

kW /nr 

Btu/in. -s 

k\V/m : 

Btu/in. : -s 

Experimental adjusted to ry* = 100 percent | 

560 

12.326 

1787.7 

3.89 

1428.00 

0.8738 

774.80 

0.4741 

413.46 

0.2530 

570 

12.645 

1834.0 

5.97 

1857.97 

1.1369 

1025.82 

.6277 

538.48 

.3295 

571 

12.488 

1811.1 

4.70 

1643.72 

1.0058 

902.76 

.5524 

470.83 

.2881 

575 

14.350 

208 1 .2 

4.65 

1815.81 

1.1 1 1 1 

1009.31 

.6176 

547.80 

.3352 

576 

14.605 

2118.2 

5.68 

2013.55 

1.2321 

1117.66 

.6839 

598.62 

.3663 

577 

14.225 

2063. 1 

4.47 

1621.66 

.9923 

942.47 

.5767 

491.58 

.3008 

580 

16.364 

2373.3 

4.27 

1856.50 

1.1360 

1070.43 

.6550 

588.98 

.3604 

TDK/MABL laminar | 

560 

12.326 

1787.7 

3.89 

536.69 

0.3284 

312.47 

0.1912 

173.23 

0.1060 

570 

12.645 

1834.0 

5.97 

745.71 

.4563 

462.00 

.2827 

255.76 

.1565 

571 

12.488 

1811.1 

4.70 

598.13 

.3660 

367.22 

.2247 

201.50 

. 1 233 

575 

14.350 

2081.2 

4.65 

630.00 

.3855 

389.77 

.2385 

216.86 

.1327 

576 

14.605 

2118.2 

5.68 

755.68 

.4624 

467.39 

.2860 

263.1 1 

.1610 

577 

14.225 

2063. 1 

4.47 

605.16 

.3703 

365.58 

.2237 

203.30 

.1244 

580 

16.364 

2373.3 

4.27 

611.37 

.3741 

384.86 

.2355 

215.23 

.1317 


Reading 

Expansion area ratio, t 


300 

388 

500 

635 

800 

975 


Heat flux to nozzle wall | 


kW/nr 

Btu/in. -s 

kW/nr 

Btu/in. -s 

kW/nr 

Btu/in. : -s 

kW/nr 

Btu/in. : -s 

kW/nr 

Btu/in. -s 

kW/nr 

Btu/in. -s 





Experimental adjusted to ty 31 

= 100 percent 





569 

297.92 

0.1823 

239.42 

0.1465 

180,91 

0.1 107 

153.13 

0.0937 

143.49 

0.0878 

122.24 

0.0748 

570 

378.33 

.2315 

297.76 

.1822 

220.79 

.1351 

170.78 

.1045 

133.84 

.0819 

111.62 

.0683 

571 

336.16 

.2057 

267.04 

.1634 

200.69 

.1228 

155.74 

.0953 

126.33 

.0773 

103.77 

.0635 

575 

384.70 

.2354 

302.50 

.1851 

225.36 

.1379 

178.30 

.1091 

140.38 

.0859 

1 18.97 

.0728 

576 

414.12 

.2534 

327.67 

.2005 

243.50 

.1490 

189.90 

.1162 

145.77 

.0892 

124.20 

.0760 

577 

353.00 

.2160 

279.95 

.1713 

209.18 

.1280 

166.20 

.1017 

136.62 

.0836 

1 18.32 

.0724 

580 

41 1.67 

.2519 

327.34 

.2003 

245.79 

. 1 504 

189.25 

.1158 

148.55 

.0909 

127.47 

.0780 

TDK/MABL laminar 

569 

117.18 

0.0717 

91.52 

0.0560 

69.62 

0.0426 

52.79 

0.0323 

40.20 

0.0246 

31.87 

0.0195 

570 

173.23 

.1060 

132.86 

.0813 

103.28 

.0632 

78.61 

.0481 

59.16 

.0362 

46.58 

.0285 

571 

137.77 

.0843 

106.88 

.0654 

80.73 

.0494 

62.59 

.0383 

46.74 

.0286 

37.59 

.0230 

575 

148.55 

.0909 

1 14.07 

.0698 

88.09 

.0539 

67.49 

.0413 

50.50 

.0309 

39.55 

.0242 

576 

178.13 

.1090 

136.13 

.0833 

106.39 

.065 1 

80.73 

.0494 

60.96 

.0373 

48.05 

.0294 

577 

138.09 

.0845 

106.72 

.0653 

83.18 

.0509 

64.23 

.0393 

47.88 

.0293 

37.91 

.0232 

580 

146.76 

' .0898 

113.42 

.0694 

87.76 

.0537 

67.00 

.0410 

50.50 

.0309 

40.20 

.0246 
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